\ 


RADC-TR-88-64 
Final  Technical  Report 

March  1988 

AD-A2G0  086 


A  THEORETICAL  STUDY  OF  THE 
EFFECTS  OF  VEGETATION  ON 
TERRAIN  SCATTERING 


Virginia  Polytechnic  Institute  and  State  University 


Dr.  Gary  S.  Brown 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griff iss  Air  Force  Base,  NY  13441-5700 


OTIC 

ECTE 

I  1 


e 


\  \  b 


tP  a  o  V.U3  ! 

V 

E 


S  8  9  22  a 


r 


'  This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Division  (PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

gADC.Xg.88.64  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 

ROBERT  J.  PAPA 
Project  Engineer 


APPROVED: 


JOHN  K.  SCHINDLER 


Acting  Director 

Directorate  of  Electromagnetics 


Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (EECE)  Hanscom  AFB  MA  01731-5000.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or 
notices  on  a  specific  document  require  that  it  be  returned. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  6f  ThB  ftAtjfe 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 
N/A 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 
N/A 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

230-11-11QF-104 -402 1511 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Virginia  Polytechnic 
Institute  &  State  University 


6c  ADDRESS  (Gty,  State,  and  ZIP  Cod*) 
Blacksburg  VA  24061 


6b.  OFFICE  SYMBOL 
(If  tpplksblt) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Rome  Air  Development  Center 


8c  AOORESS  (City.  State,  and  ZIP  Cod*) 
Hanscom  AFB  MA  01731-5000 


8b.  OFFICE  SYMBOL 
(If  applicable) 

EECE 


Form  Approved 
QMB  No.  0704-0(88 


lb.  RESTRICTIVE  MARKINGS 
N/A 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 
RADC-TR-88-64 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Rome  Ai’-  Development  Center  (EECE) 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 
Hanscom  AFB  MA  01731-5000 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

FI 9628-84 -C-001 2 


10  SOURCE  OF  FUNDING  NUMBERS 


program 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO. 

NO 

ACCESSION  NO 

61102F 

2305 

J4 

56 

11.  TITLE  (Include  Security  Cltaificadon) 

A  THEORETICAL  STUDY  OF  THE  EFFECTS  OF  VEGETATION  ON  TERRAIN  SCATTERING 


12.  PERSONAL  AUTHOR(S) 
Dr.  Gary  S.  Brown 


13a.  TYPE  OF  REPORT 

13b.  TIME  COVERED 

14.  DATE  OF  REPORT 

(Year,  Month.  Day) 

IS.  PAGE  COUNT 

Final 

FROM  Jan  85  T0  Oct  8b 

March 

1988 

68 

SUPPLEMENTARY  NOTATION 

Report  prepared  under  subcontract  from  Applied  Science  Associates,  Inc. 
Court,  Central,  SC  29630 


105  East  Windsor 


]  17.  COSATI  CODES  j 

FIELD 

GROUP 

SUB-GROUP 

20 

03 

17 

_ 02 _ 

18.  SUBJECT  TERMS  (Continue  on  revene  if  necessary  and  identify  by  block  number) 
Electromagnetic  Scattering  Multiple  Scattering 

Rough  Surface  Scattering  Asymptotic  Series 

Random  Volume  Scattering 


19.  ABSTRACT  (Continue  on  revene  if  necessary  and  identify  by  block  number) 

This  study  develops  an  analytical  model  for  the  scattering  by  a  covering  of  vegetation  on 
a  randomly  rough  conducting  surface.  The  foliage  is  analyzed  using  the  Foldy-Twersky 
formalism  for  the  coherent  or  average  field  and  a  new  distorted  wave  born  approximation 
(DWBA)  is  developed  for  the  fluctuating  field.  The  rough  surface  scattering  is  modeled 
using  an  interactive  solution  of  the  Magnetic  Field  Integral  Equation  (MFIE)  and  this 
approach  is  related  to  classical  random  surface  methods.  In  addition,  an  improvement  to 
this  interactive  approach  is  developed  based  on  a  Luneburg-Kline  series  for  the  current 
induced  on  the  rough  surface.  The  layer-on-the-surf ace  problem  is  modeled  by  a  single 
down  and  up  passage  of  the  appropriate  coherent  and  incoherent  fields  as  opposed  to  the 
exact  solution  which  comprises  an  infinite  number  of  passages.  The  key  to  the  analysis  is 
the  conversion  of  all  fields  in  the  foliage  to  either  coherent  or  incoherent  quantities. 
All  the  important  primary  and  secondary  interactions  are  preserved  in  this  model. 


20  distribution  /  availability  of  abstract 
Qunclassified/unumited  □  same  as  rpt 


□  OTIC  USERS 


22a  NAME  Of  RESPONSIBLE  INDIVIDUAL 

Robert  J.  Papa 


^nclassYftt'H 


RlTY  CLASSIFICATION 


22b  TELEPHONE  (Include  Area  Code) 

(617)  377-3585 


22c  OFFICE  SYMBOL 

RADC  (EECE) 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASS1F1CATH 


Table  of  Contents 


1.0  The  Effects  of  Vegetation  on  Terrain  Scattering 

1.1  Introduction 

1.2  Vegetation  Only 

1.2.1  The  Coherent  Field 

1.2.2  The  Fluctuating  Field 

1.2.2. 1  The  Distorted  Wave  Born  Approximation  (DWBA) 

1.2.2. 2  An  Exact  Integral  Equation  For  The  Fluctuating  Field 

1.2. 2.3  Comparison  With  Other  DWBA 

1.3  Rough  Surface  Only  -  A  Perspective 

1.4  Vegetation  On  A  Rough  Surface 

1.4.1  Summary 
References 

Appendix  A:  Integral  Equations  For  Scattering  From  Dielectric  Surfaces 
Appendix  B:  Continuation  Of  The  Coherent  Field  Into  Free  Space 
Appendix  C:  The  Conversion  of  Fields  Into  Scattering  Cross  Sections 
2.0  Asymptotic  Surface  Scattering 

2.1  Luneburg- Kline  Expansion  For  The  Surface  Current 

2.1.1  Introduction 

2.1.2  Analysis 
References 


Accession  For 

NT  IS  GRA&I 
DTIC  TAB 
Unannounced 
Justification- 


¥ 

□ 


By - 

Distribution/ 


Availability  Codes 
AvaI*L  and/or 
Special 


page 

1 

1 

2 

3 

5 

5 

7 

11 

12 

18 

30 

32 

35 

38 

41 

49 

49 

49 

51 

59 


1.0  The  Effects  of  Vegetation  on  Terrain  Scattering 
1.1  Introduction 

An  important  element  in  any  effort  to  model  the  electromagnetic  scattering  from  natural 
terrain  is  the  inclusion  of  the  effects  of  the  various  types  of  vegetative  or  foliage  cover.  That 
such  coverage  is  important  is  obvious  from  just  looking  at  terrain  and  being  able  to  readily  see 
which  surface  areas  are  covered  by  the  different  types  of  vegetation,  e.g.  crops,  forests,  scrub 
brush,  and  grass.  In  the  language  of  scattering  physics  the  phrase  "just  looking  at  terrain  and 
.  .  ."  is  equivalent  to  saying  that  the  bistatic  scattering  cross  section  of  the  vegetation  covered 
surface  in  the  visible  frequency  band  is  sufficiently  sensitive  to  the  types  of  cover  to  provide  a 
quantitative  means  for  discriminating  between  the  various  kinds  of  vegetation.  With  this  more 
technical  definition,  it  is  important  to  take  note  of  the  distinguishing  characteristics  of  the 
"measurement".  That  is,  the  transmitter/receiver  geometry  is  bistatic  and  the  transmitter  and 
receiver  are  both  operating  over  a  very  wide  band  of  frequencies.  This  latter  point  is  particularly 
important  because  such  wideband  measurements  have  no  analog  in  the  lower  frequency  ranges 
common  to  radar  or  communication,  and  so  care  should  be  exercised  in  comparing  what  is 
"seen"  with  what  much  lower  frequency  measurements  or  theory  yield.  Certainly,  one  would 
expect  that  as  the  frequency  is  continually  decreased,  there  would  come  a  point  where  the 
scattering  is  due  predominately  to  the  underlaying  surface  and  there  is  very  little  scattering  from 
the  vegetation  [1J. 

Unfortunately,  most  radars  and  communications  systems  operate  in  a  frequency  range 
that  is  somewhere  between  these  two  extremes.  Consequently,  there  is  a  need  for  models 
involving  scattering  from  (1)  the  vegetation  layer  alone,  (2)  the  underlying  surface  alone,  and 
(3)  various  degrees  of  interaction  between  the  layer  and  the  underlaying  surface.  Since  ultra 
high  spatial  resolution  scattering  is  not  a  goal  of  this  study,  the  use  of  statistical  means  to 
describe  the  scattering  should  be  appropriate.  Ideally,  this  means  that  what  is  needed  is  the 
probability  density  function  (pdf)  for  the  scattered  field  because  this  will  then  permit  a 
calculation  of  any  statistical  moment  of  the  scattered  field.  However,  the  pdf  of  the  scattered 
field  is  a  very  difficult  quantity  to  obtain,  in  general,  and  so  attention  will  be  directed  toward  a 
less  ambitious  goal,  namely,  the  modeling  of  the  first  and  second  moments  of  the  scattered  field. 

Having  established  the  goal  of  this  study,  the  next  step  is  to  develop  a  rationale  for 
achieving  it.  The  rationale  f  or  this  problem  comes  rather  naturally  from  the  importance  of  the 
two  major  contributors  to  the  scattering.  That  is,  the  first  thing  that  needs  to  be  done  is  to 
insure  that  adequate  models  exist  for  both  the  foliage  layer  and  the  underlaying  surface  alone. 
Where  deficiencies  are  noted,  they  need  to  he  corrected.  I  he  final  phase  is  to  model  the  foliage 
layer  on  the  surface  by  giving  particular  attention  to  the  interaction  of  the  two  media.  This  is 
the  real  crux  of  the  problem  because  it  is  this  phase  that  permits  a  smooth  transition  between 
the  complete  high  frequency  shadowing  of  the  surface  by  the  vegetation  to  the  low  frequency 
transparency  of  the  foliage  to  the  incident  radiation. 


1.2  Vegetation  Only 

Before  becoming  too  involved  in  the  electromagnetics  associated  with  wave  propagation 
through  and  scattering  by  vegetation,  it  is  advantageous  to  organize  vegetation  into  four  very 
broad  categories.  In  particular,  these  categories  are  trees,  brush,  crops,  and  grasses.  This  is 
done  primarily  for  discussion  purposes  because  there  is,  in  general,  little  that  is  unique  about 
each  category.  Trees  are  usually  characterized  by  their  bulk  size,  relatively  large  leaves,  and  a 
central  trunk.  Brush  is  usually  smaller  in  overall  size,  has  a  smaller  leaf,  lacks  a  well-defined 
central  trunk,  and  has  much  less  woody  structure.  However,  these  differences  become  less 
distinct  when  the  trees  are  relatively  young.  Similarly,  there  is  not  much  difference  between 
new,  emerging  grain  crops  and  grass  except  that  the  former  is  usually  planted  in  periodic  rows 
while  the  grass  tends  to  be  more  uniform  in  its  coverage.  Despite  the  biological  nonuniqueness 
of  these  categories,  they  provide  a  physically  intuitive  means  of  identifying  types  of  vegetation. 

In  estimating  the  electromagnetic  scattering  behavior  of  a  collection  of  vegetation,  it 
would  be  most  helpful  if  the  behavior  of  a  single  composite  structure  of  an  entire  entity  were 
known.  However,  even  single  trees,  bushes,  and  plants  are  difficult  to  describe  from  an 
electromagnetic  scattering  point  of  view.  Thus,  it  is  necessary  to  go  one  step  further  in 
simplicity  and  consider  the  scattering  properties  of  the  component  parts  of  the  tree,  bush,  plant, 
etc.  In  the  case  of  trees  and  bushes,  these  component  parts  comprise  leaves,  branches,  and 
trunks.  With  crops  and  grasses,  the  parts  comprise  leaves  and  stems  or  stalks.  Although  these 
component  parts  are  interconnected  by  branches,  stems,  and  twigs,  most,  subsequent  scattering 
analysis  will  ignore  these  interconnects  and  treat  the  component  parts  as  devoid  of  physical 
contact  with  any  other  parts  of  the  foliage. 

This  latter  approximation  is  clearly  a  drastic  one  and  therefore  deserves  some 
justification.  The  important  assumption  that  will  be  made  to  rationalize  this  approximation  is 
that  each  component  part  of  the  foliage  is  a  much  better  absorber  of  electromagnetic  energy 
than  it  is  a  scatterer.  More  formally,  the  approximation  is  that  the  absorption  cross  section  is 
significantly  greater  than  the  scattering  cross  section.  For  the  frequency  range  of  interest  to  this 
study,  it  turns  out  that  this  is  indeed  a  reasonably  good  assumption  [1].  Hence,  given  that  each 
component  part  of  the  vegetation  absorbs  much  more  incident  energy  than  it  scatters,  consider 
what  happens  when  a  wave  strikes  two  leaves  connected  by  a  stem  or  branch.  In  addition  to 
the  energy  scattered  from  one  leaf  to  another,  there  is  also  the  possibility  of  energy  being 
propagated  in  a  transmission  line  mode  along  the  connecting  branch  or  stem.  I  lowcvcr,  energy 
being  transported  between  leaves  via  this  mechanism  suffers  a  threefold  attenuation,  first,  it 
sulfcrs  absorption  in  the  first  leaf  and  then  in  the  connecting  branch  and  finally  in  the  second 
leaf.  Thus,  it  would  seem  reasonable  to  ignore  the  interconnection  of  the  various  parts  of  the 
vegetation  based  on  the  highly  lossy  nature  of  the  branches  and  stems,  i.e.  the  interconnecting 
medium. 


1.2.1  The  Coherent  Field 


The  total  field  inside  a  vegetative  medium  can  be  mathematically  split  into  the  sum  of 
a  mean,  average,  or  coherent  part  and  a  zero  mean  fluctuating  part  (2J,  e.g. 

Et  =  <  E{>  +  6Et  (1.1) 

where  < .  >  denotes  the  averaging  operation,  and  8Et  is  a  zero  mean  field  quantity.  There  are 
specific  physical  situations  when  analytical  models  for  the  average  field  have  been  derived.  The 
first  of  these  is  the  classic  Foldy-Twerskv  result  (3,4,5]  in  which  the  concentration  of  scattcrcrs 
is  sufficiently  small  that  multiple  scattering  between  them  may  essentially  be  ignored.  The 
Lax-Twersky  (6,7, 8, 9]  result  holds  for  higher  concentrations  of  scatterers  but  with  the  restriction 
that  they  cannot  differ  very  much  in  their  dielectric  properties  from  free  space.  Finally,  rather 
high  concentrations  of  scatterers  have  been  analyzed  recently  [10]  using  a  numerically  based 
T-matrix  approach. 

The  Lax-Twersky  theory  is  difficult  to  apply  to  the  vegetation  problem  because  the 
component  scatterers  have  such  a  large  relative  dielectric  constant.  The  numerical  T-matrix 
approach  would  be  difficult  to  apply  to  the  foliage  problem  because  of  the  diverse  shapes  of  the 
component  scatterers.  Furthermore,  neither  of  the  above  methods  is  really  necessary  because 
the  density  of  vegetation  is  relatively  low.  That  is,  the  highest  volume  fraction  occupied  by 
vegetation  will  be  less  than  5%,  and  even  this  percentage  is  relatively  unusual  (1).  More  typical 
volume  fractions  arc  1%  or  less.  With  this  low  concentration  of  scatterers  combined  with  their 
relatively  high  loss,  the  Foldv-Twerskv  theory  for  the  coherent  or  mean  field  inside  the  foliated 
medium  should  be  adequate  [11]. 

The  Foldy-Twerskv  theory  predicts  that  in  an  unbounded  volume  of  scatterers,  the 
average  field  will  be  a  plane  wave  having  the  wavenumber  kp  given  by  [12] 

*7  J  A  A 

kp  =  ko  +  *Kf  <fp(kiki)\  ni't)<Ja  (*-2) 

0 

where  k0  is  the  free  space  wavenumber.  n(<x)da  is  the  number  of  scatterers  per  unit  volume  having 
the  size  parameter  a  between  a  and  a  +  da,  and  the  total  number  of  scatterers  per  unit  volume 
is  p  where 


CO 

l>  =  f  /!<*  )da  (1.3) 

0 

A  A  , 

The  quantity  Jp(kj,kj)  is  the  pl>t  vector  component  of  the  scattering  amplitude  of  a  single 

\ 

scattcrcr  having  size  a.  It  is  computed  for  a  plane  wave  traveling  in  the  direction  k;  incident  on 


the  body  and  the  body  scattering  in  the  direction  k-v  The  polarization  of  the  incident  field  in  this 
calculation  is  taken  to  be  ep  while  the  symbol  < .  1  denotes  that  the  scattering  amplitude  may 
need  to  be  averaged  over  all  possible  orientations  of  the  scattering  body.  An  explicit  expression 
for  fp(kj,  is 

a  a  A:2 

fp(ki,ki)  =  ~^cr  ~  1 )  5  lEi(r  )'PlE‘ 1  z*P(jk0kj -T)dv  (1 .4) 

Va 

where  tr  is  the  relative  dielectric  constant  of  the  scattering  body,  £,(7)  is  the  field  inside  the  body 

A  A 

due  to  an  incident  plane  wave  having  amplitude  £',  polarization  p  ,  and  traveling  in  the  A, 
direction.  The  integration  in  (1.4)  is  over  the  volume  of  the  scatter,  Vx,  and,  alth  ugh  not 

explicitly  shown,  the  calculation  is  for  one  specific  orientation  of  the  body.  The  av  rages  in 

(1.2)  over  orientation  and  size  are  necessary  when  there  is  a  polydisperse  mixture  of  scatterers 
having  random  orientation.  As  indicated  in  (1.2),  it  is  possible  for  the  average  medium  to  be 
anisotropic,  i.e.  the  propagation  constant  is  different  in  different  directions. 

Equally  important  to  the  result  in  (1.2)  are  the  limitations  on  its  validity.  These 
limitations  were  not  particularly  well  understood  until  recently  [12]  when  it  was  shown  that  they 
could  all  be  essentially  lumped  together  in  the  following  condition; 

{kp  -  ka)  max(a)  <<  1  (1.5) 

where  max(a)  is  the  maximum  meaningful  dimension  of  the  scattering  body.  For  a  sphere, 
max(a)  is  the  diameter  while  for  a  randomly  oriented  thin  disk  it  should  be  set  equal  to  the 
diameter  of  the  disk.  The  condition  in  (1.5)  is  essentially  equivalent  to  requiring  that  the 
scattering  properties  of  a  body  in  free  space  and  in  the  random  medium  are  not  appreciably 
different.  It  should  be  noted  that  if  the  concentration  of  scatterers  increases  then  cither  the 
relative  dielectric  constant  must  approach  unity,  see  (1.4),  or  the  scatterers  must  become 
smaller,  see  (1.5),  in  order  for  (1.2)  to  remain  valid.  It  is  the  interplay  of  the  scatterer 
concentration,  dielectric  constrast,  and  size  that  determine  the  validity  of  the  Foldv-Twerskv 
result  for  the  average  or  coherent  field. 

If  there  are  a  number  of  different  distinct  types  of  scatterers  comprising  the  medium  then 
(1.2)  should  be  augmented  to  reflect  this  fact.  For  example,  if  there  are  M  different  types  of 
scatterers  then  (1.2)  should  be  written  as  follows; 

A  -)  CO  5/  A  A 

Ip  =  k~  +  4t/  1  <fJkiki)\mnm(x)d.i  (l.M 

0  '”=> 

A  A  . 

where  <  fp(kj  k-]\m  is  the  forward  scattering  amplitude  for  the  mlh  class  of  scatterers  comprising 
nm(x)dt  number  per  unit  volume  having  a  size  parameter  between  a  and  a  +  da.  This  result  is 


particularly  relevant  to  the  vegetation  problem  because  there  arc  indeed  a  number  of  distinctly 
different  types  of  scatterers  (even  when  the  actual  foliage  components  are  replaced  by  canonical 
shapes).  Table  I  shows  the  major  vegetation  categories  that  have  been  postulated  along  with 
the  actual  component  parts  comprising  these  and  the  canonical  approximations  to  these  parts. 
The  canonical  shapes  are  the  ones  that  would  be  used  in  calculating  scattering  amplitudes  such 
as  required  by  (1.6). 

Clearly,  there  is  a  great  deal  of  approximation  involved  in  replacing  exact  shapes  by 
canonical  forms.  However,  it  must  be  remembered  that  vegetation  is  a  very  complicated 
physical  environment  and  it  is  essential  that  it  be  simplified  as  much  as  possible  without 
distorting  the  important  details.  This  was  the  approach  used  in  an  earlier  attempt  to  go  beyond 
pure  empirical  modeling  [1]  and  it  appeared  to  work  quite  well  for  wave  propagation  through 
trees.  Since  this  work  first  appeared,  there  has  been  a  great  deal  of  literature  on  computing  the 
scattering  properties  of  various  canonical  shapes  and  it  seems  that  there  are  a  number  of 
approximations  that  can  be  used  [13,14,15].  These  results  are  most  welcome  because  calculating 
the  scattering  amplitudes  of  even  the  limited  number  of  essential  canonical  shapes  is  not  a  trivial 
matter. 

1.2.2  The  Fluctuating  Field 

Equation  (1.1)  indicates  that  the  total  field  inside  the  random  medium  comprises  a 
coherent  or  average  field  and  a  zero  mean  fluctuating  part.  The  coherent  field  attenuates 
exponentially  due  to  absorption  and  scattering  by  the  objects  in  the  medium;  the  attenuation 
rate  is  determined  by  the  imaginary  part  of  (1.4).  The  fluctuating  field  gives  rise  to  the 
incoherent  scattered  power  and  is  therefore  of  prime  importance  to  this  study.  The  fluctuating 
field  is,  in  general,  more  difficult  to  determine  than  the  coherent  field;  however,  reasonable 
success  has  been  obtained  using  the  so-called  distorted  wave  Born  approximation  (DWBA)  [16] 
at  least  for  the  frequency  range  of  interest  to  this  study.  However,  as  the  DWBA  is  now  used 
[17],  there  appears  to  be  a  form  of  double  accounting  for  the  effects  of  the  average  medium. 
To  show  this,  the  DWBA  will  first  be  reviewed,  an  exact  integral  equation  for  the  fluctuating 
field  will  be  derived,  and  this  result  will  be  compared  to  the  DWBA  that  is  presently  in  use.  The 
differences  will  be  noted  and  discussed. 

1.2. 2.1  The  Distorted  Wave  Born  Approximation  (DWBA) 

Taylor  [16]  gives  an  excellent  discussion  of  the  DWBA  as  it  applies  to  particle  scattering 
by  potentials.  1  louevcr,  the  same  reasoning  can  be  applied  to  the  vegetation  scattering  problem 
and  the  rationale  for  the  approximation  centers  around  the  form  of  ( 1 . 1 ).  What  (1.1)  docs  is  to 
split  the  unknown  field  into  the  sum  of  a  largely  dominant  term.  <  If  >,  and  a  small 
perturbation  <5 Et  .  1  he  dominant  term  <  If  >  is  known  to  a  reasonably  good  approximation; 
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BASIC  VEGETATION 
CATEGORY 


Tree 


: 


Brush 


Crops 


Grasses 

! 


ACTUAL  COMPONENT  CANONICAL  COMPONENT 
PARTS  SHAPES 


leaves/needles 

branches 

trunks 


disks/cylinders 

cylinders  with  large  range  of  radii 
cylinders  with  small  range  of  radii 


leaves 

branches 


disks 

cylinders 


leaves 

stalks,  stems 
product 


disks  on  thin  elongated  sheets 
cylinders  with  relatively  small  range 
of  radii 
variable 


leaves 

stems 


short  dipoles 
cylinders 


TABLE  I 


A  categorization  of  basic  vegetation  types  by  actual  and  canonical 

component  parts. 


hence,  6Et  can  be  determined  fairly  accurately  by  a  rather  crude  approximation.  The  key  to  the 
success  of  the  DWBA  is  that  most  of  Et  is  already  known  ( <  Et  > )  and  it  is  only  a  slight 
improvement  (<5 E()  that  is  sought.  Of  course,  this  also  shows  one  of  the  limitations  of  the 
DWBA  as  applied  to  the  vegetation  problem.  If  <  Et  >  is  to  be  the  dominant  term  in  Et ,  the 
DWBA  cannot  be  expected  to  give  an  accurate  estimate  for  SEt  very  deep  into  the  medium. 
This  is  because  <  Et>  will  be  very  small  for  deep  penetration  into  the  medium  and  it  is  then  that 
SEt  will  become  the  dominant  field  quantity.  This  situation  clearly  violates  the  assumptions  in 
the  DWBA. 

In  the  next  section,  an  exact  integral  equation  of  the  second  kind  for  the  fluctuating  field 
will  be  developed  and  attention  will  be  directed  toward  the  Born  term  in  this  equation.  Using 
this  term  only  comprises  a  DWBA  because  it  assumes  that  <  Et>  is  known  and  dominant.  In 
the  subsequent  section,  this  result  will  be  compared  to  a  heuristic  based  DWBA.  It  will  be 
shown  that  the  heuristic  approach  leads  to  what  appears  to  be  a  form  of  double  accounting  for 
the  average  properties  of  the  medium.  In  this  regard,  the  heuristic  approach  is  not  in  agreement 
with  the  Born  term  from  the  exact  expression  for  the  fluctuating  field. 


1.2. 2. 2  An  Exact  Integral  Equation  for  the  Fluctuating  Field 

The  total  field  at  the  point  r  ,  Et(r  ),  due  to  an  incident  field,  £j(r  ),  and  the  fields  scattered 
from  objects  located  about  the  point  r ,  Es(r),  is  given  by 


Et{r)  =  Ej(r)  +  Es(r) 


(1-7) 


The  scattered  field  due  to  the  presence  of  the  N  objects  can  be  written  as  follows  [18]; 


Es{r  )  =  LK-LEt(rn  +  rQ) 


(1.8) 


where  L  is  the  three  dimensional  integral  over  all  space,  i.e.,  L  —  f f  fdr0  ,  and  the  dyadic 
operator  K  £  is  given  by 


*> 


,v  _  =  _  _ 

V  (e  -  1  )Sn(r0\  QnW{r  -  r„  -  r()) 
n=  1 


(1.9) 


In  ( 1.9)  cr^  is  the  relative  dielectric  constant  of  the  n ^  scattercr.  Sn(r()\  i1n)  is  the  support  of  the 
n scattercr  whose  centroid  is  located  by  the  position  vector  rn  so  that 


■s>,r 


1  r  inside  I  n 

0  rQ  not  inside  l  'n 


(1.10) 


and  V  is  the  volume  of  the  n^1  scattercr.  The  dyadic  T  is  given  by 


/ 


(1.11) 


r  =  -p.v\j  +  k0h0i^S{r  -  7n  -  70)  +  is(7  -  7n  -  70)iikl 

where  P.V.  denotes  the  principle  value*,  7  is  the  unit  dyad,  S( . )  is  the  three-dimensional  delta 
distribution,  and 

g(r  ~  rn~  ro)  =  exp (~jk0\7  -  7n  -  rQ\)l4n\7  -  7n-70\  (1.12) 

The  variable  Sdn  in  the  argument  of  the  support  function  symbolizes  the  dependence  on  the 
orientation  of  the  n1^  scatterer.  Equation  (1.9)  is  general  in  that  it  allows  scattcrer-to-scatterer 
variation  in  the  dielectric  constant  (er^)  ,  the  size  and  shape  (K^),  the  orientation  or  alignment 
(ii, j),  and  the  location  ( rn )  .  It  should  be  noted  that  calculation  of  the  scattered  field  requires 
knowledge  of  the  total  field  inside  each  scatterer.  Thus,  substituting  (1.8)  in  (1.7)  and  then 
taking  the  point  r  inside  each  scatterer  yields  N  coupled  integral  equations  for  the  total  field 
inside  each  of  the  N  objects.  Once  these  equations  are  solved,  the  scattered  field  at  any  point 
in  space  can  be  determined  by  straightforward  integrations.  For  the  random  problem,  a 
somewhat  difTerent  methodology  will  be  developed. 

Substituting  (1.8)  into  (1.7)  yields 

Et  =  Et  +  LKZE{  (1.13) 

The  total  field  at  any  point  in  space  can  be  written  as  the  sum  of  a  conditional  average,  <  ZTf], 
and  a  zero-mean  fluctuating  part,  6 Et,  so  that  (1.13)  becomes 

<  E[\  -f-  <5£f  —  Ej  4-  LK^Lj  (1.14) 

This  decomposition  on  the  left  side  is  a  bit  different  from  the  one  usually  used  and  results  from 
the  need  to  accommodate  the  following  situation  unique  to  discrete  object  scattering.  The 
conditional  average  <  £  ]  is  an  average  of  the  total  field  with  the  point  of  observation  of  the 
field  held  constant,  i.e.  the  point  of  observation  is  the  conditioned  variable.  Under  most 
circumstances,  the  point  of  observation  is  a  constant  and  so  the  conditional  average  equals  the 
total  average.  This  would  be  the  case  on  the  left  hand  side  of  (1.14),  e.g.  <  £;]  =  <  Et>. 
Furthermore,  the  conditionally  averaged  total  field  is  the  one  that  is  usually  called  just  the 
’mean''  nr  "average"  or  "coherent"  J'teld.  However,  it  is  the  riylit  side  of  (1.14)  that  gives  rise  to 
a  need  for  the  conditionally  averaged  total  field.  That  is.  the  random  quantities  in  the 
expression  for  Av  ,  1 1.9).  are  usually  the  scuttercrs'  positions  (rn)  si/cs.  (I)(),  and  orientations 

2-\s  Jt  appears  in  (1.11),  the  P.V.  excludes  a  small  spherical  volume  centered  at 


S 


(Q„)  .  The  field  quantity  on  the  right  hand  side  of  (1.14)  has  an  explicit  dependence  on  the 
coordinates  rn  +  r0\  rn  locates  the  centroid  of  the  n1^1  scatterer  while  r0  ranges  over  its  support 
Vn.  Thus,  the  scatterer  dependent  terms  on  the  right  side  of  (1.14)  may  be  written  as  follows; 

L%Et  =  lFx(T  -  7n  -  70)  Et  (rn  +  70)  (1.15) 

When  Et(rn  +  r0)  is  decomposed  into  a  sum  of  a  conditional  average  and  a  zero  mean 
fluctuating  term,  there  results 

Et(7n  +  70)  =  <  Et(70  +  rn)]  4-  6Et(7n  +70)  (1.16) 

If  both  sides  of  (1.16)  are  averaged,  it  is  done  by  first  performing  a  conditional  average  over  all 
random  variables  except  rn  .  Such  a  conditional  average  of  SEt  yields  zero  because  this  is  the 
way  6Et  was  formed,  i.e.  <  <5£,]  =  0.  The  total  average  of  Et  thus  amounts  to  averaging 

<  Et(rn  +  rQ)1  over  the  random  variable  7n  . 

The  point  of  the  above  discussion  is  to  illustrate  the  need  for  using  a  conditional  average 
rather  than  a  total  average  in  decomposing  Et .  If  a  total  average  had  been  used,  there  would 
have  been  no  way  of  determining  what  to  do  about  the  rn  dependence  on  the  right  side  of  (1.14). 
When  using  a  conditional  average  it  becomes  clear  what  to  do  both  when  the  point  of  field 
observation  is  random  and  nonrandom.  In  the  latter  case,  the  conditional  average  is  equal  to 
the  full  average. 

Substituting  (1.16)  in  the  rhs  of  (1.14)  yields 

<  E[\  +  <5£j  —  £)  4-  iA’v<  E[j  4-  LKydEf  (1.17) 

Averaging  this  equation  gives 

<  =  Ej  +  L  <  K £  <  E^\  >  4 -  L  <  K >  (1.18) 

and  subtracting  (1.18)  from  (1.17)  leads  to  the  following  integral  equation  for  6Er; 

<)E[  =  L(  1  —  P)K  v  <  ZTj]  +  L(  1  —  P)k  v<)/Tf  (1-19) 


where  P  is  the  averaging  operator,  i.e. 


At  this  stage,  (1.19)  is  an  exact  integral  equation  for  the  fluctuating  component  of  the  total  field. 
However,  the  fluctuating  part  of  the  total  field  is  equal  to  the  fluctuating  part  of  the  scattered 
field.  This  can  be  proved  by  writing  (1.7)  as 

<£,]  +  6Et  =  £)  +  <£j  4-  6ES  (1.22) 

averaging  this  equation  and  subtracting  the  result  from  (1.22)  to  give 

6E,  =  SEs  (1.23) 


Thus,  (1.19)  becomes 


6ES  =  L(  1  -  P)KZ  <  £f]  +  L(l  -  P)KzSEs  (1.24) 

which  is  now  an  integral  equation  for  the  fluctuating  part  of  the  scattered  field.  The  Born  term 
in  (1.24)  depends  on  the  average  of  Et. 

The  average  or  mean  part  of  the  Born  term  in  (1.24)  is  thus  given  by 

LPKZ  <  £,]  =  fK z(r  -Tn-  70)p(rn)  <  E,(fn  +  70)]d7n  (1.25) 

n—  1,2, ... ,  ;V 


which  is  a  well  defined  quantity. 

If  all  of  the  scattercrs  are  confined*  to  the  volume  V,  equation  (1.24)  relates  the 
fluctuating  scattered  field  in  this  volume  to  the  fluctuating  scattered  field  at  any  point  in  space. 
If  the  point  of  observation  of  SEs  on  the  left-hand  side  of  (1.24)  is  taken  inside  V  then  (1.24) 
becomes  an  integral  equation  for  6ES  (inside  V).  Solving  this  integral  equation  inside  V 
subsequently  allows  the  determination  of  6ES  outside  V  through  the  use  of  (1.24).  The  Born 
term  in  (1.24)  depends  on  the  fluctuating  part  of  the  product  of  the  propagator  and  the  averaged 
total  field.  It  is  interesting  to  note  that  there  are  no  average  or  effective  medium  quantities 
appearing  in  (1.24)  other  than  <  Et ]  .  The  Born  term  should  be  dominant  whenever  <  £fJ  docs 
not  depart  too  much  from  free  space  propagation.  When  <  £r]  violates  this  condition,  it  then 
becomes  necessary  to  account  for  the  second  term  in  (1.24).  For  foliage  and  vegetation  in  the 
frequency  range  of  interest  to  this  study,  the  Born  term  in  (1.24)  should  be  adequate  for  most 
all  cases,  i.e. 


This  is  equivalent  to  saying  that  the  position  vector  n~  1,  2 . V  locating  the  centroid 

of  the  nt/l  scattcrcr  must  iie  within  the  volume  V  -  where  depends  on  where  the  centroid 
of  the  nt/l  scattcrer  is  located  (within  the  scattcrcr). 
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(1.26) 


5ES  as  U 1  -  /5)/CI<  £,] 

Even  though  this  result  comes  from  the  Born  term  of  an  exact  integral  equation,  it  can  be 
thought  of  as  a  distorted  wave  Born  approximation  (DWBA)  because  of  the  assumption  of  the 
dominant  nature  of  the  coherent  or  mean  field  [16].  Clearly,  the  approximation  of  (1.26)  merits 
further  study  primarily  because  it  comes  from  an  exact  integral  equation. 


1. 2.2.3  Comparison  With  Other  DWBA 

There  are  other  forms  of  the  DWBA  that  have  been  applied  to  foliage.  Lang  [17,13]  has 
developed  a  heuristic  form  of  the  DWBA  in  which  the  scatterers  are  immersed  in  a  medium 
having  the  same  effective  dielectric  constant  as  the  mean  or  average  field,  e.g. 


<tr>p  =  kjlkj,  (1.27) 

where  k£  is  given  by  (1.2).  The  scatterers  are  then  assumed  to  be  illuminated  by  <  Et]  and  to 
scatter  independently  of  each  other.  Thus,  any  multiple  scattering  effects  included  in  this  model 
are  contained  entirely  in  the  use  of  (1.27)  as  the  background  medium  and  <  Et ]  as  the  incident 
field.  If  the  analysis  technique  derived  in  the  previous  section  is  applied  to  this  scattering 
situation,  the  resulting  Born  term  corresponding  to  (1.26)  is 


-  /V  _  _ 

6ES  =  m-nz.  <  Et~] 
n- 1 


(1.28) 


where 


~  kp(Crn  *)^/i(ro;  kp(r  rn  r<f  ( 1 .29) 

is  the  propagator  for  a  scatterer  of  relative  dielectric  constant  er^  buried  in  a  medium  having  a 
wavenumber  equal  to  kp.  Thus,  /Cj  is  the  average  propagator  for  the  n scatterer. 

In  comparing  the  heuristic  result  of  (1.29)  with  the  exact  Born  term  in  (1.26),  there  is 
only  one  essential  difference  and  that  is  the  use  of  an  effective  background  medium  rather  than 
free  space.  It  may  be  possible  that  the  heuristic  approach  comprises  a  partial  summation  of 
some  of  the  contributions  from  the  second  term  on  the  rhs  of  (1.24):  however,  this  could  be 
dillicult  to  prove.  On  the  other  hand,  the  early  and  pioneering  work  of  Twcrsky  [4,5.7]  appears 
to  argue  against  the  use  of  an  effective  medium  propagator  and  in  favor  of  the  free  space 
propagator.  In  fact,  the  use  of  an  effective  medium  propagator  seems  closely  akin  to  double 
accounting  for  the  mean  field  effects.  There  was  not  enough  time  to  resolve  this  discrepancy 
between  the  exact  and  heuristic  approaches.  However,  it  is  recommended  that  this  point  be 
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investigated  further  since  the  two  results  lead  to  different  answers  and  the  difference  may  be 
appreciable  for  relatively  thick  foliage  layers.  A  potentially  fruitful  approach  to  reconciling  the 
two  results  is  to  try  some  form  of  selective  partial  summation  of  the  second  term  on  the  rhs  of 
(1.24).  That  is,  solving  (1.24)  via  iteration  yields 

6ES  =  Z  0(1  -  P)KI)m [L(l  -  P)KZ  <  £f]l  (1.30) 

m=0  L 

or 

6ES  =  m-/>)f2  I  {L(l-P)KI)m<E[]  (1.31) 

m— 0 

and  it  may  be  possible  to  manipulate  this  latter  series  in  such  a  fashion  as  to  derive  an 
approximate  Born  term  which  uses  from  (1.29)  rather  than  K  v  . 

1.3  Rough  Surface  Only  -  A  Perspective 

For  the  problem  of  wave  propagation  through  and  scattering  by  discrete  random  media, 
there  are  a  rather  limited  number  of  analytical  techniques  for  dealing  with  this  problem  [19]. 
Conversely,  when  analyzing  scattering  by  randomly  rough  surfaces,  there  seem  to  be  a  great  deal 
of  diverse  methods  capable  of  producing  results  (20).  The  key  word  here  is  "seem"  because,  in 
fact,  many  of  the  methods  lead  to  essentially  the  same  results;  the  major  differences  arc  with  the 
starling  point  and  the  degree  of  approximation  necessary  to  achieve  the  result.  In  view  of  this 
observation,  it  would  be  highly  desirable  if  a  fundamental  approach  could  be  found  which,  when 
subjected  to  a  hierarchy  of  approximations,  could  be  shown  to  lead  to  a  similar  hierarchy  of 
approximate  results.  The  primary  advantage  resulting  from  finding  such  a  technique  is  that,  as 
the  level  of  the  approximation  is  reduced,  there  is  a  guaranteed  increase  in  the  accuracy  of  the 
resulting  solution.  Consequently,  there  is  a  clear  reason  for  pursuing  improvements  in  this 
technique  because  they  will  result  in  an  improved  solution. 

The  technique  that  appears  to  hold  the  greatest  potential  for  leading  to  a  rigorous  solution 
to  the  rough  surface  scattering  problem  is  the  Magnetic  Field  Internal  Equation  (MFIE)  for  the 
current  induced  on  a  perfectly  conducting  surface.  The  infinite  conductivity  limitation,  for  the 
frequencies  of  interest  to  this  study,  is  significant  only  near  grazing  incidence.  Appendix  A 
considers  the  case  when  the  surface  or  interface  is  a  dielectric.  Essentially  what  happens  is  that 
there  is  both  an  electric  and  a  magnetic  current  induced  on  the  surface.  These  currents  satisfy 
coupled  integral  equations  of  the  second  kind;  the  kernels  of  the  equations  arc  more  complicated 
than  with  the  MFIE,  but  the  basic  equations  arc  the  same  as  the  MITE.  These  complications 
arc  not  unexpected  as  it  is  known  that  the  dielectric  surface  may  give  rise  to  effects  which  are 
not  observed  with  the  conducting  surface.  For  example,  a  Brewster  angle  phenomenon  with  the 
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coherent  scattered  field  should  be  possible  along  with  a  reentry  effect  in  which  the  incident  field 
enters  one  side  of  a  bump  on  the  surface  and  exits  on  the  other  side.  As  noted  above,  the  effects 
are  only  significant  near  grazing  incidence  and  they  are  very  dependent  on  the  imaginary  part 
of  the  dielectric  constant  of  the  surface.  The  main  point  of  this  discussion  is  that  if  techniques 
can  be  developed  for  the  MF1E  which  are  not  strongly  dependent  on  the  particular  form  of  the 
kernel  then  it  should  be  possible  to  translate  the  bulk  of  the  methodology  to  the  dielectric 
surface  integral  equations.  As  noted  above,  this  topic  is  discussed  more  fully  in  Appendix  A. 

Our  analysis  of  the  MF1E  as  it  applies  to  rough  surface  scattering  is  a  continuing  effort. 
For  example.  Section  2.1  presents  a  technique  which  extends  the  geometrical  optics  solution  of 
the  MFIE  into  the  nonzero  wavelength  regime.  However,  it  seems  appropriate  at  this  point  to 
provide  some  perspective  to  previous  analytical  solutions  of  the  MFIE;  that  is,  to  show  where 
they  come  from  and  how  they  fit  in  a  well-ordered  hierarchy  of  approximations.  As  noted  above, 
this  is  essential  if  we  are  to  establish  the  fact  that  the  MFIE  leads  to  successively  more  accurate 
solutions  as  the  input  approximations  are  similarly  improved. 

The  form  of  the  MFIE  that  is  most  familiar  is  the  following: 

Jjr)  =  (O  x  H‘(7)  +  2/T (7)  X  /  Js(r0)  x  V0<7(  |7  -7Q\  )dsQ  ( 1 .32) 

where  Js  is  the  surface  current  density,  n  is  the  unit  normal  to  the  surface,  //'  is  the  incident 
magnetic  field,  and  G  is  the  free  space  Green's  function 

G(i r  -  rQ |)  =  exp  ( -jkQ\7  -  r0|)/4rr|7  -  7*  |  (1.33) 

Equation  (1.32)  can  be  rewritten  in  the  following  form 

7(7)  =  27(7)  x  /?(7 )  +  /  [  -  7(7)  x  V 0G  x7(70)]  d7to  (1.34) 

where 

7  =  jsk  i  +  c2  +  9 1/2 
7  =  -  ;rc  -  c/  +  2 

ds0  =  dr,f\  +  c:  +  9 
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and  dr,  =  dx0dy0.  The  Cx  =  dCldxa/idt^y  =  <?C Idy  are  the  x  and  y  components  of' the  surface 
slope.  The  right  hand  side  can  be  further  manipulated  by  standard  vector  identities  to  yield 

7(7)  =  2N  x  If‘  +  /  {  -27flG  [7(7) .  +  2[7(7)  •  7 0G]  ]7  (r0)<Ir,^  (i.3f 


where 


70g[7(7).  )7  =  v0g[/V(7).7(70)] 

In  operator  notation  (1.35)  can  be  written  as  follows: 


J  -  Ji  +  LK.J 


where 


7'  =  2 7  x  /?’ 


(1.38a) 


L  =  /  (  ) 


(1.386) 


K  =  {-27 0G  [7(7).  +  2[/V(7).V0G]} 


(1.38c) 


Integral  equations  of  the  second  kind  such  as  (1.37)  can  be  formally  solved  via  iteration  to  yield 


J  =  Z.(LK-)nJj 
n= 0 


This  scries  can  be  thought  of  as  a  series  of  iterates,  i.c. 


7  =  £  7n) 

n= 0 


where 


7m)  =  {LK.)m7j  (1.41) 

and  the  /croth  order  iterate  =  ./,■  is  called  the  Born  term.  The  scattered  electric  field 
resulting  from  the  current  Js  is  given,  in  the  far  field  approximation,  as  follows: 


£r(  R0)  =  jk(), 1()  G(  R„)ks  xktxf  ./t(  r  )  exp  ( ik t .  r  )ds 
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where  the  vector  R0  points  from  the  origin  of  the  coordinate  system  on  the  mean  surface  in  the 

A  —“A 

direction  ks  to  a  distance  R0  away.  Also,  the  vector  ks  is  defined  by  ks  =  k0ky  The  integral  in 
(1.42)  can  be  written  as 


Es(Ro)  =  jko1oG(Ro)ks  x  ks  x  fJ  (r)  exP  Hksp  exp  Uks(  •~rt)drt  ( 1  •4-5) 

where  rj0  is  the  impedance  of  free  space  and 

ks  =  kS{  +  ksz 

The  integral  in  (1.43)  is  recognized  to  be  the  2-dimensional  Fourier  transform  of  the  current  J 
weighted  by  the  exponential  factor  exp  (jks  £),  i.e. 


W  =  Jko^oG(Ro>ks  x  ks  xF2^-J  exP  VksS)] 


where 


F2[  •  ]  =  /  •  exp  (jkS{  .  rt)drt. 
Expanding  the  double  curl  operation  leads  to 

W  =  jk0n0G(Ro)F&J L  exP Ukspy 


where 


J±  =  J  ~  ks^ks  ‘ J ) 

— *  .  ,  A 

is  the  component  of  J  which  is  perpendicular  to  kr  Assuming  that  the  iterative  series  in  (1.40) 
converges  so  that  the  series  can  be  integrated  term  by  term  yields 

ES(R0)  =Jk0r,0G(R0 )  I  ?w)  (1.44) 

m= 0 


where 


-p[m) 


5"> 


exp  (iksp 


( 1 .45) 


The  above  development  shows  that  there  are  two  important  elements  to  the  determination 
of  the  scattered  field.  First,  the  current  iterates  must  be  computed  from  ( 1.41 )  and,  second,  the 
corresponding  scattered  field  iterates  must  be  calculated  from  (1.45).  Of  course,  what  must  be 
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done  and  what  is  actually  done  arc  two  entirely  different  situations  and  this  difference  leads  to 
various  degrees  of  approximation.  In  the  following  material,  an  ordering  of  most  of  the  classical 
approximations  will  be  developed  based  upon  their  effect  on  the  current  and  scattered  field 
iterates.  The  purpose  of  this  ordering  is  to  put  these  classical  approximations  into  a  perspective 
which  clearly  indicates  where  we  have  been  and  were  we  need  to  proceed. 

The  first  level  of  approximation  is  low  frequency  in  nature.  It  starts  with  exp (jk.s£) 

as  an  approximate  form  of  the  current  and  further  simplifies  this  to 

exp (jksp  *  2Hj  x /i;{  1  +  j{ks^  -  exp  (-  jk (1.46) 

where  the  slopes  in  7(0)  have  been  ignored  and  the  surface  roughness  is  assumed  to  be  so  small 
that 


exp [j{ks  -  kj  )C]  %  1  +  j{ks  -  kj  ){  (1-47) 

z  ^  z 

The  other  factors  in  (1.46)  come  from  the  form  of  the  incident  magnetic  field,  i.e. 

Hl  =  HJn  exp  (  -jki£  -  jkl(.7() 

Integrating  (1.46)  exactly  to  form  leads  to  a  scalar  Bragg  result  in  the  backscatter  direction 
for  the  incoherent  scattered  power  (21].  This  result  follows  from  the  fact  that  7^)  in  the 
backscatter  direction  is  polarized  in  the  same  direction  as  the  incident  electric  field,  and  that  it 
depends  on  £  in  a  linear  fashion. 

The  next  level  of  approximation  is  to  use  an  exact  form  for  7(0)  but  evaluate  £j0)  by 
stationary  phase  techniques  so  that  the  result  is  a  high  frequency  limiting  form.  This  solution 
is  frequently  called  specular  point  scattering  [22].  It  should  be  noted  that  whereas  the  scalar 
Bragg  solution  results  from  an  exact  integration  of  an  approximate  form  for  7(0),  the  specular 
point  result  comes  from  an  approximate  integration  of  the  exact  7(0).  The  specular  point  result 
assumes  that  ( ks  -  kjJZ  >  >  1,  the  surface  slopes  are  relatively  small,  and  the  surface  curvatures 

*■  i. 

are  also  small  (large  radii  of  curvature).  The  first  approximation  permits  the  use  of  stationary 
phase;  the  second  eliminates  multiple  scattering  on  the  surface;  the  last  avoids  sharp  edge 
diffraction  effects. 

A  final  approximation  involving  only  7(())  applies  to  composite  surfaces.  Such  surfaces 
have  a  range  of  scales  which  satisfy  both  the  large  scale  approximation  common  to  specular 
point  scattering  and  the  small  scale  approximation  to  the  scalar  Bragg  result  |2.'|.  If  the  surface 
can  be  approximated  by  a  composite  surface,  then  the  integration  to  give  l,\' 1 1  can  be 
accomplished  essentially  exactly.  This  leads  to  the  composite  surface  scattering  result  in  which 
the  specular  point  term  dominates  around  the  specular  direction  and  the  scalar  Bragg  term  is 
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important  at  other  angles.  In  addition,  the  scalar  Bragg  result  is  tilted  by  the  large  scale  surface 
slopes.  It  is  important  to  remember  however  that  the  scalar  Bragg  result  shows  no  polarization 
sensitivity. 

A  summary  of  what  the  zeroth  order  iteration  of  the  MFIE  leads  to  under  certain 
approximations  is  given  in  Figure  1.1.  It  should  be  noted  that  there  is  no  discussion  of  when 
is  a  valid  approximation  for  the  total  current.  This  question  can  be  answered  only  by  using 
improved  estimates  of  the  current  to  compute  the  scattered  field.  What  is  done  above  and 
summarized  in  Figure  1.1  is  to  show  the  consequences  of  (a)  approximating  and  doing  the 
integration  to  get  £$°)  exactly,  (b)  reversing  this  situation,  and  (c)  using  a  combination  of  these 
manipulations. 

The  exact  same  procedure  may  be  repeated  when  the  next  iteration  of  the  current  is 
included.  That  is,  the  current  J  is  next  approximated  by 

7  =  70)  4-  (1.48) 

where 


7°)  - 


and 


71}  = 


LK.Jt 


It  has  recently  been  shown  [24]  that  a  low  frequency  approximation  ofyl°)  4-  /*)  along  with 
the  assumption  of  small  surface  slopes  leads  to  a  result  which  when  integrated  exactly  to  yield 
4-  E gives  rise  to  the  vector  Bragg  result.  The  vector  Bragg  result  differs  from  the  scalar 
solution  obtained  with  jw)  only  in  that  for  backscatter,  the  vector  solution  shows  a  polarization 
sensitivity  whereas  the  scalar  does  not.  This  illustrates  the  inadequacy  of  the  /**)  iterate  (or 
Born  term)  in  so  far  as  the  scattered  field  in  the  o(T-specular  direction  is  concerned. 

The  next  level  of  approximation  involves  evaluating  all  integrations  in  the  high  frequency 
limit.  Thus,  when  LK  •  Ji  is  so  evaluated  it  leads  to  -  Jj  in  the  shadowed  parts  of  the  surface, 

0  on  the  illuminated  parts  of  the  surface,  and  the  possibility  of  first  order  multiple  scattering 
[25|.  That  is.  first  order  multiple  scattering  means  the  refection  of  an  incidence  ray  from  one 
point  of  the  surface  to  another  point.  Generally,  this  first  order  ray  optic  multiple  scattering  is 
ignored  because  it  makes  the  scattered  field  integration  very  hard  to  do.  That  current  which 
remains  is  the  shadowed  current,  i.c.  +  A'f  If  the  integrations  needed  to  compute 
+  /:{ 1  ^  are  evaluated  asymptotically  (as  !<0—oo),  what  results  is  the  shadowed  specular  point 
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solution.  Comparing  this  solution  with  the  corresponding  limiting  behavior  of  £^°)  alone  shows 
that  the  inclusion  of/1)  in  this  limit  gives  rise  to  shadowing. 

Finally,  if  the  surfaces  under  study  can  be  split  into  both  large  scale  structure  (Cs)  and  small 
scale  undulations  (£5)  so  that  high  frequency  techniques  can  be  used  to  predict  the  scattering 
from  the  large  structure  and  low  frequency  methods  can  be  used  on  the  small  scale  structure 
then  the  surface  can  be  treated  as  a  composite  structure,  i.e.  C  =  £<f  +  Cs  •  One  can 
subsequently  show,  though  not  easily  (26j,  that  if  those  parts  of/0)  +  Z1)  which  depend  upon 
are  treated  by  high  frequency  asymptotics  and  those  parts  which  depend  on  are  treated 
by  low  frequency  approximations  then  the  result  is  shadowed  specular  point  scattering  plus 
tilted  vector  Bragg  diffraction.  The  tilting  of  the  Bragg  scattering  is  by  the  large  scale  surface 
slopes.  It  is  interesting  to  note  that  the  sum  of  only  /°)  and  Z1)  lead  to  one  of  the  best 
scattering  models  around,  namely,  the  composite  scattering  model. 

The  results  of  including  both  Z°)  and  Z1)  are  summarized  in  Figure  1.2.  In  all  categories 
there  is  a  marked  improvement  over  using  Z°)  only.  However,  this  is  particularly  so  for  the  low 
frequency  results  where  there  is  a  polarization  sensitivity  not  accounted  for  with  Z°)  only. 

The  purpose  of  this  section  has  been  to  put  the  various  levels  of  rough  surface  scattering 
approximations  into  perspective.  This  has  been  done  by  relating  the  approximations  to  either 
the  first  two  iterate  solutions  of  the  MFIE  or  the  subsequent  scattering  integral  calculations. 

It  has  been  shown  how  this  methodology  can  lead  very  naturally  to  a  hierarchy  of 
approximations.  Most  important  of  all,  however,  is  the  fact  that  the  MFIE  can  generate 
improved  solutions  by  simply  increasing  the  order  of  the  iteration.  This  establishes  the  MFIE 
as  a  useful  means  for  studying  rough  surface  scattering. 

Finally,  all  of  the  above  approximate  solutions  are  available  for  the  surface  scattering  part 
of  the  vegetation  layer  problem.  Which  one  should  be  used  is  dictated  by  the  frequency  of 
interest  and  the  surface  roughness  statistics. 

1.4  Vegetation  on  a  Rough  Surface 

Section  1.2  and  1.3  detail  the  techniques  that  can  be  used  to  estimate  the  scattering  from 
an  isolated  patch  of  vegetation  and  bare  ground,  respectively.  This  section  will  develop  and 
discuss  a  model  for  scattering  from  the  combination  of  a  layer  of  vegetation  on  a  rough  surface. 

Some  of  the  earliest  work  on  this  problem  proposed  using  the  simple  addition  of  the 
scattering  cross  sections  of  the  vegetation  layer  and  the  rough  surface.  However,  it  was  quickly 
realized  that  some  accounting  for  the  inllucncc  of  the  vegetation  layer  on  the  rough  surface 
scattering  cross  section  was  essential.  Subsequent  attempts  to  account  for  the  vegetation  simply 
attenuated  the  surface  cross  section  by  the  loss  suHcred  by  the  coherent  power  propagating 
down  through  the  foliage  and  then  back  up  to  the  foliage  air  interface.  L'sc  of  this 
approximation  led  to  some  qualitative  agreement  with  scattering  measurements,  but  the  model 
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Bragg 


Figure  1.1  The  description  of  the  scattered  field  resulting  from  the  zeroth  order  iteration  of  the 
MFIE  current,  coupled  with  various  approximations. 


Vector  Bragg 
(No  Tilting) 


Incident  Shadowed 
Specular  Point 


Shadowed  Specular  Point 
+  Tilted  Vector  Bragg 


Figure  1.2.  The  description  of  the  scattered  field  resulting  from  the  zeroth  and  first  order 
iteration  of  the  MFIE  current,  coupled  with  various  approximations. 
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also  raised  questions  as  to  why  it  worked  as  weil  as  it  appeared  to  and  when  it  could  be  expected 
to  fail.  Very  little  of  this  early  work  was  ever  published  and  most  of  this  information  was 
presented  during  informal  meetings  and  discussions. 

More  recently,  attempts  were  made  to  rigorously  model  certain  subsets  of  this  general 
problem  (see  [  1 3 j  for  a  comprehensive  list  of  relevant  references).  Unfortunately,  most  of  these 
attempts  chose  to  model  the  foliage  as  a  continuous  variation  in  dielectric  constant  and  to  thus 
use  continuous  random  media  theories  as  the  analytical  tools.  While  this  approach  did  lead  to 
some  insight  into  the  scattering  problem,  it  also  raised  questions  because  of  the  need  for 
"effective"  or  "equivalent"  parameters  in  the  continuous  representation  of  the  discrete  random 
medium.  Lang  and  Sidhu  [13]  overcame  the  limitations  of  these  earlier  models  by  using  the 
Foldy-Twersky  theory  along  with  a  version  of  the  DWBA  to  model  the  scattering  from  a  foliage 
layer  on  a  flat  earth.  This  was  a  significant  advance  because  it  clearly  showed  the  interaction 
between  the  foliage  components  and  the  reflection  from  the  Hat  surface.  Unfortunately,  it  is 
not  obvious  how  Lang  and  Sidhu's  method  of  analysis  could  be  extended  to  an  arbitrarily 
roughened  surface. 

The  model  that  will  be  developed  here  is  very  simple  in  concept  but,  as  one  might  expect, 
complicated  in  detail.  First,  the  field  incident  on  the  foliage  is  converted,  inside  the  foliage,  to 
a  mean  or  average  field  and  a  zero  mean  fluctuating  field.  The  mean  field  is  determined  by  the 
Foldy-Twersky  theory  of  Section  1.2.1  while  the  fluctuating  field  is  based  on  the  distorted  wave 
Born  approximation  of  Section  1.2.2. 1.  These  two  fields  then  act  as  incident  fields  on  the  rough 
surface  which  is  approximated  as  perfectly  conducting.  The  electric  currents  excited  by  these 
fields  are  computed  using  the  Magnetic  Field  Integral  Equation  (MF1E)  as  discussed  in  Section 
1.3.  The  next  step  is  to  let  these  currents  radiate  back  up  through  the  foliage.  The  fields 
radiated  by  the  surface  currents  are  then  acted  on  by  the  foliage  to  give  rise  to  another  set  of 
mean  and  fluctuating  fields  in  the  foliage  [27],  These  fields  can  be  continued  into  free  space  via 
standard  techniques  such  as  volume  and  surface  distributions  of  current  or  plane  wave  spectral 
methods.  If  the  process  is  truncated  at  this  point  then  what  has  been  accounted  for  is  one  single 
downward  and  upward  passage  of  the  fields.  If  the  process  were  to  be  continued,  the  interaction 
of  the  upward  going  fields  with  the  foliage  should  be  allowed  to  interact  with  the  rough  surface 
which,  in  turn,  would  scatter  back  up  through  the  foliage.  If  the  problem  were  to  be  analyzed 
exactly,  this  process  should  be  continued  an  infinite  number  of  times.  However,  there  are  very 
practical  reasons  for  not  continuing  this  process  beyond  the  first  downward  and  then  upward 
pass.  The  most  obvious  of  these  is  the  computational  complexity  associated  with  each 
down-and-up  field  iteration.  Furthermore,  each  iteration  will  also  generate  the  need  lor  higher 
order  (multipoint)  vegetation  and  surface  statistics  and  these  arc  simply  not  known.  Finally,  if 
there  is  an  indication  that  more  down-and-up  field  iterations  arc  needed  then  perhaps  there  arc 
simpler  techniques  available  to  deal  with  the  entire  problem.  The  single  down-and-up  field 


20 


technique  assumes  that  either  the  foliage  or  the  surface  is  the  dominant  scattercr  and  that 
interactions  between  the  two  natural  scatterers  is  small  relative  to  these  two  dominant  effects. 
To  a  certain  extent,  this  approximation  can  be  checked  and  this  will  be  discussed  later  in  the 
section. 

The  foliage  and  the  rough  ground  represent  statistically  different  random  processes.  If 
there  were  some  degree  of  correlation  between  the  two,  it  would  be  necessary  to  use  conditional 
averaging  in  forming  the  first  two  moments  of  the  scattered  field.  That  is,  the  foliage  would  be 
averaged  holding  the  surface  fixed  and  then  the  surface  would  be  averaged.  This  process 
requires  the  use  of  conditional  probability  density  functions.  However,  there  is  no  reason  to 
expect  to  a  first  order  at  least  that  there  should  be  any  correlation  between  the  foliage  and  the 
surface.  In  fact  we  will  go  further  and  assume  that  the  two  processes  are  statistically 
independent;  this  simplifies  the  actual  mathematical  operations  because  the  joint  density 
function  is  just  a  product  of  marginal  densities.  Thus,  the  averaging  over  either  the  foliage  or 
the  surface  can  be  done  independent  of  the  remaining  random  process.  However,  it  should  be 
remembered  that  if  there  is  some  correlation  between  the  foliage  and  the  surface  then  it  can  be 
accommodated  in  the  model  through  the  use  of  conditional  probability  density  functions.  This 
generality  may  be  useful  in  dealing  with  very  high  resolution  scattering  because  such  a  situation 
may  emphasize  the  slight  correlations  between  certain  classes  of  foliage  and  terrain. 

The  terrain  is  assumed  to  be  a  perfectly  conducting,  randomly  rough  interface,  having  a 
zero  mean  about  the  z  =  0  plane  and  homogeneous  but  arbitrary  statistics.  The  foliage  or 
vegetation  occupies  the  space  immediately  above  the  surface  and  up  to  an  average  height  of 
z  -  h.  It  should  be  remembered,  however,  that  this  average  height  of  the  foliage  is  a  description 
which  must  be  somewhat  carefully  interpreted.  For  example,  any  coherent  scattering  which 
results  from  this  average  planar  height  is  clearly  open  to  question  since  the  foliage  interface,  for 
coherent  calculations,  should  be  modeled  as  a  rough  interface.  This  roughness  will  certainly  give 
rise  to  a  very  highly  attenuated  coherent  field  which  can  thus  be  ignored.  This  is  not  to  say  that 
there  will  never  be  a  coherent  field  scattered  from  a  discrete  collection  of  scatterers;  there  clearly 
will  be  such  a  field  when  the  volume  fraction  of  objects  approaches  unity.  However,  for  a 
sparsely  populated  medium,  the  existence  of  a  coherent  scattered  field  would  appear  to  be  a 
model  anomaly  rather  than  a  physical  reality. 

The  field  incident  on  the  foliage  layer  is  assumed  to  have  the  following  plane  wave  form: 

f:)  =  [  flvcv  +  /:/,(■/,]  exp  (  -//',•  •  r  )  (  l.4‘») 

where 

kj  =  k{  x  +  kjy  +  kjj  (1.50) 

x  v  r 
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and 


k-t  =  —  kQ  sin  0 1  cos  4>j 
kj \  =  —  k0  sin  Oj  sin  4>j 
kj^  =  -kQ  cos  6 i 

The  unit  vectors  and  denote  vertical  and  horizontal  polarization  directions; 


and 


A  ,  A  ,  A 

c/,  =  —  sin  <p  jX  +  cos  <p  jy 


(1.51a) 

(1.516) 

(1.51c) 


(1.52a) 


=  —cos  Oj  cos  <j>jX  —  cos  Oj  sin  4>jy  +  sin  Of  (1.526) 

The  angles  are  standard  polar  (Oj)  and  azimuthal  (<£,)  angles  and  they  are  explicitly  defined  in 
Figure  1.3.  As  discussed  above,  we  ignore  the  possibility  of  any  coherent  rellection  from  the 
foliage  layer.  Thus,  according  to  Section  1.2,  the  in  'idem  field  is  converted  to  the  sum  of  a 

— t.  — 

coherent  field,  <  Ej-d>  ,  and  a  fluctuating  field,  6Efd,  inside  the  foliage.  The  average  field  in  the 
foliage  is  given  by 


<  Efd>  =  E/V  exp  (  -jk^r  )  +  £/,?/,  exp  (  •  r  )  (1.53) 

where  the  directions  of  k ^  and  kj  are  the  same  as  kj ,  i.e. 

kiy  =  ki!ko  4  =  Wko 

but  the  complex  amplitudes  are  solutions  of  (1.2)  with  p=  v  or  h,  e.g. 

2  2  A  A 

kv  =  k0  +  AnJ  <fv(kj,kj)]fTlnm(x)dat.  (1.54) 

h  0  h 

and  </v{kj,kj)]m  is  the  scattering  amplitude  of  the  mUl  component  of  the  foliage  averaged  over 
all  possible  orientation  angles.  As  discussed  in  Section  (1.2),  the  integration  in  (1.54)  accounts 
lor  a  size  distribution  of  the  mlk  foliage  components.  The  fact  that  the  average  field  may  have 
different  wavenumbers  or  propagation  constants  for  horizontal  (h)  and  vertical  ( v)  polarizations 
has  been  observed  in  some  measurements  (1|;  however,  the  effect  seems  to  disappear  above 
about  800  \|  f  [/. 
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The  fluctuating  field  SEj-j  is  a  solution  of  the  integral  equation  in  (1.19);  however,  based 
on  the  discussion  presented  in  Section  1.2.2  it  should  be  possible  to  use  only  the  Born  term  in 
(1.19)  and  this  is  equivalent  to  the  distorted  wave  Born  approximation  (DWBA),  i.c. 

SEfdxL(\-P)Kz<Efd]  (1.55) 

where  L  is  the  integral  over  all  space,  P  is  the  averaging  operator,  K  £  is  given  by  (1.9),  and 
<  Ejtf]  is  essentially  the  average  field  evaluated  at  ~r0.  The  reader  is  reminded  that  a  complete 
discussion  of  the  terms  in  (1.55)  and  the  approximations  under  which  it  is  valid  are  given  in 
Section  1.2.2.  It  should  be  emphasized  that  all  of  the  quantities  are  known  and  so  (1.55)  is 
computable.  There  are  two  additional  calculations  that  need  to  be  done  with  the  fluctuating 
field  in  (1.55).  First,  it  is  a  primary  source  of  incoherent  power  scattered  back  up  into  free 
space,  i.e.  <  |<5£/y|2>.  In  fact,  in  the  absence  of  the  terrain  surface  or  for  a  very  thick 
vegetation  layer  this  is  the  only  source  of  scattered  power.  The  second  calculation  recognizes 
that  SEj-j  and  <  Ej-J]  are  both  incident  on  the  terrain  surface  and  they  induce  currents  on  the 
surface  which  reradiate  back  into  the  vegetative  medium. 

Within  the  limitation  discussed  in  Section  1.3,  the  terrain  will  be  approximated  as  a 
perfectly  conducting  interface.  Consequently,  the  Magnetic  Field  Integral  Equation  (MFIE) 
describes  how  the  induced  current,  J  ,  behaves  on  the  surface,  i.e. 

7  =  2N  x  Hfd  +  2/V  x/7  x  V0Gd7tg  (1.56) 

The  incident  magnetic  field  Hj-  may  be  written  as 

11  fd  =  +  811  fd 

and 

<  "fdl  =  -  — l—*  x  <  £fd\  <l-58fl) 

<3 llrj  =  — 7— — 7  x  SEfj  (1.586) 

'  ja}H0 

follow  from  Maxwell's  equations.  Since  <  Ej-j |  is  a  plane  wave  traveling  in  the  ^-direction 
through  a  medium  having  an  average  relative  dielectric  constant  <  rT>p  =  (p  =  v  or  h), 

( 1.58a)  becomes 

1/2 

-  (  <  V  >P)  A  - 

<  Ufd\  =  - - - ki  x  <  I[fd\  0-59) 
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where  r\Q  is  the  impedance  of  free  space  (Vm0/e0  )•  The  choice  of  vertical  or  horizontal 
polarization  in  (1.59)  determines  both  the  polarization  of  <  Ej-f]  and  <  cr>p.  If  near  zone  and 
Fresnel  zone  fields  are  ignored  in  the  calculation  of  6Ej-d  and,  subsequently,  611 yd  then  (l.5Sb) 
can  be  written  as  follows; 


»«fit  -  -£*  X  &fi  (1.60) 

A  — *■  A 

where  k  specifies  the  direction  of  propogation  of  6Efd.  It  is  not  necessary  that  k  point  in  only 
one  direction  as,  for  example  it  would  in  the  case  of  a  plane  wave.  The  average  medium  relative 
dielectric  constant  does  not  appear  in  (1.60)  because  according  to  the  DWBA  used  in  (1.55),  the 
fluctuating  field  propagates  in  free  space,  this  point  was  discussed  in  detail  in  Section  1.2. 2. 3. 
It  should  be  noted  that  all  the  quantities  in  the  Born  term  in  (1.56)  are  known. 

Substituting  (1.57)  into  (1.56)  yields 

7  =  2/V  x  <  Hjf]  +  2 Nx6Hj-d  +  2 Nx  f  7x7 QGdrt  (1.61) 

The  total  current  J  is  split  as  follows; 

J  =  Jaf  +  J6f  ( i  .62) 

where 

J, af  =  2N  x  <  Ufd  +  2 iVx  /  7af  x  V 0G^1q  ( 1 .63) 

and 

=  2N  x6Hj-d  +  2Nxf7df  xV 0Gd7,o  (1.64) 

From  (1.63),  ^-is  the  current  induced  on  the  surface  by  the  downward  propagating  average 
foliage  field  <  Hj-f]  while,  from  (1.64),  J^j- is  the  current  induced  on  the  surface  by  the  downward 
traveling  fluctuating  foliage  field  &Hfd-  These  equations  can  be  solved  by  iteration,  such  as 
discussed  in  Section  1.3,  to  yield 


Jaf=  I  (/A/)m(2.Vx  <//,/]) 

m=() 

( 1  •65') 

7sf=  I  (L.\/)W(2.7 x <5 ///•/) 

m=() 

(1.66) 

where  L  is  two-dimensional  integral  over  the  z  =  0  plane  and 


(1.67) 


.1/  =  -  270Gf  -  'gov  (7) .  +  2[,v  (7 ) .  v0G(7  -  g]} 

The  quantities  in  (1.67)  are  defined  in  Section  1.3. 

The  currents  in  (1.65)  and  (1.66)  are  now  taken  to  reradiate  upwards  into  free  space.  The 
calculation  of  the  resulting  fields  can  be  accomplished  via  either  the  exact  formula, 

Esu=  -  fj-1  xl  xf  J  (fQ)G{R -7Q)drt  (1.68) 

Ko  u 

or  the  far  field  approximation, 

-*•  A  A  A_,_* 

Esu  K  jko^oG(R)ksyksxJ  J  (ro)  exP  (Jkoks  •  r0>drt0 -  ( 1  -69) 

or  a  plane  wave  spectral  representation  developed  in  Appendix  B.  The  easiest  of  these  to  deal 
with  is  the  far-field  approximation  of  (1.69)  but,  of  course,  this  does  not  make  it  correct. 
However,  there  do  not  appear  to  be  any  near-field  effects  which  would  be  augmented  or 
magnified  by  the  foliage  layer  so  the  far-field  approximation  will  be  used.  The  total  field 
scattered  by  the  surface  and  up  into  free  space  may  be  written  as 

Esu  =  Easu  +  Edsu  (1-70) 


where 


Easu  ■  jk0r!oG(R)ks^fJaj(7o)  exp  {jkQks .  7a)drt(_ 


(1.71) 


and 


A  A  — *  A  ^  k 

E6su  =  jkot1oG(R')ksxks*f  J6j(ro)  exP  Ukoks-ro')drt0  l1-72) 

The  quantity  Easu  is  the  electric  field  scattered  up  into  free  space  by  the  rough  surface  when 
illuminated  by  the  average  foliage  field  <  Hff\.  E$su  is  the  electric  field  scattered  up  into  free 
space  by  the  rough  surface  when  illuminated  by  the  fluctuating  foliage  field  o/lfj. 

The  next  step  in  the  process  is  to  let  the  fields  in  (1.71)  and  ( 1.72)  be  incident  upon  the 
foliage  layer  from  below.  The  reason  lor  doing  this  is  that  this  problem  can  be  treated  just  like 
the  first  part,  i.e.  the  downward  passage  of  the  free  space  incident  field  through  the  foliage. 
That  is.  the  total  field  in  the  foliage  due  to  Eaw  and  E-)su  incident  from  below  can  be  written  as 
follows; 


26 


(1.73) 


p-sf  -  (^asu  +  pdsu>  + 

This  equation  is  the  same  as  (1.13)  but  with  the  incident  field  replaced  by  (E^  +  E$su).  We 
now  develop  the  standard  method  of  smoothing  approach  as  a  means  to  find  the  average  and 
fluctuating  upward  going  fields  in  the  foliage.  All  averages  are  over  the  random  quantities  in 
the  foliage  and  not  the  surface;  averaging  over  an  ensemble  of  surfaces  will  be  delayed  until  later 
in  the  analysis.  First,  Esj- in  (1.73)  is  rewritten  in  terms  of  its  average  and  fluctuating  parts,  i.e. 


<  Esj\  +  <5 Esf=  E^  +  E6su  +  LK^k  Esj\  +  LKl5Esj~  (1.74) 

Averaging  this  equation  yields 

<  fy]  =  Easu  +  LP *T  <  Esj\  +  LPKz6Esf  (1.75) 

and  subtracting  (1.75)  from  (1.74)  gives 

5Esf  =  ESsu  +  L(l  -  P)Kz<Esf]  +  HI-  P)Ki»Esf  0-76) 

Assuming  for  the  moment  that  <  Esj]  is  known  allows  the  integral  equation  in  (1.76)  to  be 
solved  by  iteration  to  yield 

Spsf  =  p6su  +  £  (Ml  -  PWzf  C  <  Esfl  +  5  J  O'77) 

n=  1 


Comparing  the  integral  equation  in  (1.76)  with  the  corresponding  equation  for  the  first 
downward  pass  through  the  foliage,  i.e.  eqn.  (1.24),  it  is  noted  that  (1.76)  has  the  additional 
term  E^su.  This,  of  course,  is  the  field  scattered  upward  by  the  surface  when  illuminated  by  the 
downward  propagating  fluctuating  field.  Thus,  to  the  first  pass  scattered  fluctuating  field  given 
by  (1.55)  must  be  added  the  scattering  of  5L ^  from  the  rough  surface,  i.e.  E$su.  It  should  be 
noted  that  both  SEj-j  and  pSsu  propagate  in  free  space  and  not  in  the  average  medium. 

In  order  to  find  the  total  fluctuating  field  due  to  the  surface  scattered  fields,  <  Esj]  must 
be  determined.  To  do  this,  (1.77)  is  substituted  in  (1.75)  for  <5£5yand,  after  regrouping  terms, 
gives 

<  ESA  =  EdW  +  LPKy  I  {/„(  1  -  P)Ky}%su  +  LPKy  I  {/.(l  -  P)Ky},1<eA  ( 1 -7S) 

J  “i=0  *=<)  “ 

which  is  the  desired  integral  equation  for  <  Esj].  In  order  to  understand  the  meaning  and 
implications  of  this  equation,  it  is  beneficial  to  write  the  corresponding  equation  for  the  mean 


1 
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field  propagating  downward  through  the  foliage.  This  topic  was  covered  in  Section  1.2.1,  but 
the  integral  equation  for  the  mean  field  was  not  given.  It  is  as  follows; 


<Ef]  =  Ei  +  LPKzj:  {U  1  -  P)Kz)n<Ej\ 


(1.79) 


where  Ei  is  the  field  incident  in  free  space.  Comparing  (1.78)  and  (1.79)  shows  that  the  only 
difference  is  the  Born  term  and  this  is  as  it  should  be  because  the  eigenvalues  (kp)  for  the  mean 
field  should  be  independent  of  downward  or  upward  wave  travel.  What  this  comparison  means 
is  that  the  average  field  propagating  upward  from  the  surface  and  through  the  foliage  is  given, 
within  the  framework  of  the  Foldy-Twersky  approximation,  by  the  Born  term  in  (1.78)  with  k0 
in  (1.71)  and  (1.72)  replaced  by  kp  (the  wavenumber  for  the  average  field).  That  is, 


<  v3  -  eL  +  s  iw  -  Vi/4U 

*  n=0  K 

where 

P-asu  =  jkp*ipOk(R)ksxksxfJaj(ro)  exp  ( jk[jls.70)drl \q 
ESsu  =  jkpi1pGk(R)KxksxIJ6^o>  exP  (/'Vj' 


(1.80) 


(1.81) 


(1.82) 


As  a  reminder,  <  Esj]  is  the  upward  propagating  field  in  the  foliage  after  having  been  scattered 
by  the  surface  and  averaged  over  all  possible  foliage  configurations.  A  further  simplification 
of  (1.80)  comprises  taking  only  the  n  =  0  term  in  the  scries,  i.e. 


<Es/]  *E%SU  +  LPKZfElsu  (1.83) 

This  result  has  a  number  of  interesting  features  which  deserve  comment.  First,  the  field  E%su  is 
the  obvious  contribution.  It  is  the  field  produced  by  the  surface  scattering  of  the  downward 
propagating  field  back  up  into  the  average  medium  characterized  by  the  wavenumbers  kp,  p  = 
v  or  h.  This  is  why  all  of  the  quantities,  except  Jaf,  in  (1.81)  depend  on  kp  rather  than  kiy,  i.e. 
E^u  is  scattered  into  the  average  medium.  The  current  Jaf,  as  it  appears  in  (1.81),  should  also 
be  computed  using  kp  rather  than  k0.  however,  such  a  replacement  is  necessary  onlv  when  there 
is  significant  multiple  scattering  on  the  surface  because  it  aticmiatcs  the  multiple  scattering.  In 
view  of  the  fact  that  we  have  no  surface  scattering  theories  to  deal  with  strong  multiple 
scattering,  there  is  no  point  to  replacing  k()  by  kp  in  computing  the  surface  current. 

The  last  terms  on  the  rhs  of  (1.83)  is  deceiving  because  it  appears  that  we  are  obtaining  a 
coherent  result  front  an  incoherent  field.  The  field  lih ...  is  indeed  a  zero  mean  field  as  can  be 
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seen  from  (1.82)  and  (1.66);  however,  the  average  of  the  product  of  A'^  and  E^su  is  not  zero. 
The  cause  of  this  is  a  highly  cooperative  effect  in  which  the  fluctuating  field  scattered  by  a  given 
object  is  rescattered  by  the  surface  back  up  to  the  same  object  which,  in  turn,  scatters  the  field 
a  third  time  and  back  into  free  space.  The  very  unique  aspect  of  this  scattering  arrangement  is 
that  the  total  length  of  the  scattering  path  is  a  constant  regardless  of  where,  in  the  random 
foliage  layer,  the  scatterer  is  located.  This  is  not  true  if  there  is  only  a  single  scattering  from  an 
object,  such  as  causes  6 Ej-j,  or  if  there  is  scattering  between  different  objects.  This  particular 
source  of  coherent  scattering  is  not  expected  to  be  large  compared  to  £ ^  or  when  the 
randomness  of  the  surface  is  averaged  over;  however,  there  are  instances  where  it  might  be 
important.  To  the  author's  knowledge,  this  effect  has  not  been  previously  noted. 

Before  leaving  the  upward  propagating  average  field  in  the  foliage,  as  given  by  (1.83),  we 
must  get  this  field  into  the  free  space  above  the  foliage  layer.  If  the  frequency  is  below  about 
400  MHz  then  there  is  a  potential  for  a  lateral  wave  existing  at  the  interface  because  <  £Ty]  is 
propagating  into  a  less  dense  medium  and,  in  this  case,  the  techniques  developed  by  Tamir  |28) 
should  be  used  to  complete  the  analysis.  However,  since  the  frequency  of  interest  to  this  study 
is  considerably  above  400  MHz,  the  lateral  wave  can  probably  be  ignored  and  the  fields  can  be 
continued  up  into  free  space  through  the  use  of  a  simple  Huygen  s  source  on  the  z  =  h  plane  or 
the  application  of  the  plane  wave  spectral  approach  developed  in  Appendix  B. 

Finally,  we  turn  our  attention  to  the  upward  traveling  fluctuating  field  as  given  by  (1.77). 
Truncating  the  series  in  (1.77)  with  the  n  =  1  term  yields 


It  is  interesting  to  compare  (1.84)  with  (1.55)  which  is  the  downward  propagating  fluctuating 
field  resulting  form  the  first  pass  through  the  foliage.  The  field  E$su  is  the  result  of  the  surface 
scattering  6 £^  into  the  upward  direction.  The  reaction  of  the  foliage  components  to  this  field 
is  given  by  the  term 


£(1  -  P)KzE6su 

and  the  remaining  term  in  (1.84)  is  due  to  the  usual  interaction  of  the  average  field  with  the 
scattcrcrs  such  as  in  (1.55).  The  important  point  to  note  is  that  the  effect  of  an  incident 
fluctuating  field  is  twofold;  it  contributes  by  itself  and  by  interacting  with  the  scattcrcrs. 
Substituting  (1.83)  into  (1.84)  yields 

^sf*  +  *■<»-  +  F-Osu)  (1^) 
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The  interesting  point  about  this  result  is  that  6Esj-  includes  fields  which  propagate  in  free  space 
and  fields  which  propagate  in  the  average  medium  up  to  a  particle  and  then  they  are  scattered 
into  free  space. 

1.4.1  Summary 

This  section  has  been  concerned  with  the  development  of  a  model  for  the  field  scattered 
by  a  foliage  or  vegetation  layer  on  a  rough  surface.  Emphasis  has  been  placed  on  the  foliage 
aspects  of  the  problem  because  they  are  the  most  difficult  to  properly  account  for  and  the 
assumed  statistical  independence  of  the  foliage  and  the  surface  permits  dealing  with  the  surface 
averaging  at  any  point  in  the  problem.  The  analysis  is  based  on  the  way  in  which  a  foliage  layer 
converts  an  incident  field  into  a  coherent  or  average  field  and  a  zero  mean  fluctuating  field. 
We  use  this  particular  decomposition  to  track  the  field  down  through  the  foliage  and  onto  the 
rough  surface.  In  the  processes  of  doing  this,  we  determine  the  field  scattered  by  the  foliage 
components  back  up  into  free  space.  This  is  the  field  that  would  exist  if  the  foliage  depth  were 
infinite. 

The  integral  equation  for  the  current  induced  on  the  rough  surrace  is  solved  by  iteration 
and  the  fields  rescattered  back  up  into  the  foliage  layer  is  expressed  in  terms  of  a  standard 
diffraction  integral.  The  source  fields  for  the  up-going  problem  included  both  coherent  and 
incoherent  terms  in  contrast  to  the  down-going  problem  which  had  only  a  deterministic  field 
as  the  source.  These  fields  led  to  some  additional  interactions  not  present  in  the  down-going 
problem.  The  net  results  were  double  pass  average  and  fluctuating  fields.  The  fluctuating  fields 
could  be  continued  directly  into  the  space  above  the  foliage  layer  because  in  the  DWBA  they 
propagate  in  free  space  even  in  the  foliage  layer.  The  average  or  coherent  up-going  fields 
require  a  Huygen  s  source  or  a  plane  wave  spectral  decomposition  at  the  foliage-air  interface 
because  they  propagate  in  an  average  medium  while  inside  the  foliage  layer. 

In  view  of  the  volume  of  equations  developed  in  this  section,  it  is  useful  to  summarize  the 
important  ones.  The  results  have  been  left  in  terms  of  field  quantities  and  no  surface  averaging 
has  been  undertaken. 

The  fluctuating  field  scattered  by  the  foliage  components  on  the  downward  passage  of  the 
incident  field  is  given  by 


di:Jd  =  U  1  -  /’)AV  <  (I  SM 

where  L  is  the  integral  over  all  space,  P  is  the  averaging  operator.  Av  is  given  by  (1.9),  and 
<  P.j-j]  is  given  by  ( 1.53)  with  r  replaced  by  rQ.  This  field,  <5/y^.  is  valid  in  all  space  above  the 
rough  surface. 
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The  fluctuating  field  scattered  by  the  through-the-foliage  passage  of  the  surface  scattered 
field  is  given  by 

4/  =  +  u- 1  -  $su  +  w  ( 1  -87) 

where  is  given  by  (1.72),  E\ m  is  given  by  (1.81),  is  given  by  (1.9)  with  k0  replaced  by 
kp  as  given  in  (1.54),  and  E^su  is  given  by  (1.82).  This  field,  SEsj-,  is  also  valid  in  all  space  above 
the  rough  surface. 

Finally,  the  average  or  coherent  field  propagating  up  through  the  foliage  from  the  rough 
surface  is  given  by 

=  (1.88) 

This  field  is  valid  only  in  the  foliage  layer  and  a  Huygen  s  source  or  a  spectral  decomposition 
is  needed  to  continue  these  fields  into  the  space  above  the  foliage  (see  Appendix  B). 

The  above  relationships  are  the  primary  results  of  this  study.  It  is  recommended  that  they 
be  used  to  develop  numerical  scattering  results  which  can  be  checked  against  measured  radar 
data. 
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Appendix  A 

Integral  Equations  For  Scattering  From  Dielectric  Surfaces 

When  trying  to  determine  the  electric  current  induced  on  a  perfectly  conducting  body  or 
surface,  there  are  three  integral  equation  formulations  that  are  familiar.  The  first  is  an  integral 
equation  of  the  second  kind  called  the  Magnetic  Field  Integral  Equation  (MF1E).  The  second 
is  an  integral  equation  of  the  first  kind  called  the  Electric  Field  Integral  Equation  (EFIE). 
These  equations  are  based  on  the  discontinuity  or  continuity  of  the  tangential  components  of 
the  stated  fields  on  the  conducting  boundary.  A  third  equation  is  of  the  first  kind  and  it  results 
from  the  requirement  of  zero  total  field  inside  the  conducting  medium.  Of  these  equations  or 
formulations,  we  have  chosen  the  MFIE  because  a  formal  solution  by  iteration  can  always  be 
developed.  In  fact  this  was  done  in  Section  1.3  and  the  significance  of  the  first  two  iterates  was 
discussed.  Convergence  of  the  iterative  series  is  always  a  potential  problem  but  there  is  also  the 
possible  solution  by  partial  summation  methods.  This  is  why  we  favor  equations  of  the  second 
kind  over  integral  equations  of  the  first  kind. 

When  the  surface  becomes  an  imperfect  dielectric  or  a  lossy  conductor,  there  are  many 
possible  ways  to  describe  the  field  scattering  problem.  A  good  discussion  of  the  possible 
formulations  is  given  by  Jones.*  In  addition  to  volume  and  surface  integration  formulations, 
he  discusses  the  use  of  integral  equations  developed  from  the  Leontovich  impedance  boundary 
condition.  This  latter  formulation  is  approximate  but  very  attractive  because  one  can  deal  with 
either  an  electric  or  a  magnetic  current  in  a  single  integral  equation.  There  are  two  difficulties 
with  the  technique.  First,  the  technique  is  approximate**  and  its  ability  to  describe  scattering 
from  arbitrarily  roughened  surfaces  is  as  yet  unknown.  Second,  the  actual  resulting  equation 
is  an  integro-difFerential  equation  in  that  it  involves  the  unknown  current  and  its  derivative. 
For  these  reasons,  integral  equations  based  upon  the  impedance  boundary  condition  do  not 
seem  to  be  suitable  to  the  random  rough  surface  problem. 

In  searching  for  integral  equations  suitable  for  describing  the  induction  of  currents  on  the 
dielectric  interface,  wc  were  guided  by  one  primary  condition.  This  condition  was  that  the 
equations  be  sufficiently  close  to  the  MFIE  that  we  could  use  MFIE-based  methods  to  deal 
with  the  dielectric  interface.  That  is.  wc  did  not  want  to  deal  with  the  dielectric  interface  as  a 
fundamentally  new  and  uilferent  problem.  A  number  of  formulations  were  investigated  and  the 

*  Jones,  D.S.,  Methods  It.  Electromagnetic  (fair  Propagation.  Oxford  Press,  pp.  SS7,  1979. 

**  Wang,  D.S.,  "Limits  and  Validity  of  the  Impedance  Boundary  Condition  on  Penetrable 
Surfaces",  IEEE  Trans.  Antennas  <£  Propag.,  AP-JJld).  pp.  453-4.57,  1987. 
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one  developed  by  Muller***  was  selected  as  the  most  suitable.  In  this  approach  Muller  derives 
coupled  integral  equations  of  the  second  kind  for  the  electric  current  Js  and  the  magnetic  current 
Ks.  The  resulting  equations  are  as  follows; 

7sf)  =  (7~Ti~)^(r  }  +  (^  +  1)  ~  ^oG^dso 

-j-Jl  +  k~\Kf0) .  70][70(70  -  V 0G])ds0  (. A 1) 

W)  =  (7^7)<(0  +  )x/Ks{70)x[V0G0  -  cr70G]ds0 

+  ~  +  k-\7f0) .  Vj[V^0  -  7  0G]}ds0  (AT) 

where  tr  and  n r  are  the  relative  permittivity  and  permeability  of  the  lower  medium, 
=  fv0lco  is  characteristic  impedance  of  free  space,  k0  =  2rr/A0,  k  -  and 

7S  =  nxfi1  K‘s  =  7‘xn  (A3) 

G0  =  exp  (  -  jkQ  |7  -  7q  |  )/4>r  |7  -  7q  |  (/14) 

G  =  exp  (  — yVc | r  -7*  |)/4t|7 -7*|  (.45) 

The  operator  70  is  the  conventional  three  dimensional  gradient  operator  evaluated  on  the 
surface. 

Of  particular  note  with  these  equations  is  that  they  are  second  kind  integral  equations  and 
they  involve  only  the  unknown  currents  and  not  their  derivatives.  It  is  only  necessary  to  solve 
(Al)  for  Js x  and  Js^  and  (A2)  for  Ks x  and  Ks •  because  the  currents  must  be  tangential  to  the 
surface.  That  is,  n  •  Js  =  0  and  n  .  Ks  -  0  so  that 

Jsz  =  W,  -  l>JSy  (-<6) 

and 

K*z  =  ;*Ksx  * 

***  Muller.  D..  Foundations  of  Mathematical  Theory  of  Electromagnetic  Springer- Vcrlag, 

pp.  l')6‘). 


where  £v  and  (y  are  the  x  and  y  components  of  the  surface  slope.  Thus,  (Al)  and  (A2)  can  be 
rewritten  as  a  matrix  integral  equation  of  the  second  kind  for  K$x  and  KSy.  While  the 

increase  in  dimensionality  relative  to  the  conducting  interface  case  is  not  insignificant,  it  does 
not  cause  any  fundamentally  new  problems.  The  major  increase  in  difficulty  comes  from  the 
need  to  deal  with  the  second  derivative  of  the  Green's  functions  at  the  source  point.  However, 
this  is  also  is  not  a  fundamentally  new  problem. 

The  point  of  the  above  discussion  is  to  show  that  the  penetrable  rough  interface  problem 
can  be  dealt  with  using  a  set  of  equations  that  are  very  similar  to  the  MFIE  describing  the 
conducting  interface.  This  means  that  if  techniques  can  be  developed  to  deal  with  the  MFIE 
which  are  not  tailored  to  the  very  special  form  of  the  kernel  of  the  integral  equation  then  they 
can  also  be  applied  to  the  penetrable  interface  also.  Of  course,  special  kernel  dependent 
techniques  can  also  be  translated  from  one  problem  to  the  other  as  long  as  both  kernels  have 
the  same  property.  In  summary,  it  is  the  similarity  of  the  MFIE  with  (Al)  and  (A2)  which 
justify  our  emphasis  on  developing  analytical  techniques  for  solving  the  MFIE. 


Appendix  B 

Continuation  of  the  Coherent  Field  Into  Free  Space 

Section  1.4.1  summarizes  the  various  contributions  to  the  scattered  field  in  the  region  of 

space  above  the  foliage  layer.  These  contributions  comprise  coherent  and  fluctuating  fields  and 

the  latter  present  no  difficulty  because  they  are  easily  continued  into  free  space  (see  the 

discussion  in  Section  1.4.1).  However,  the  coherent  or  average  fields  propagate,  when  in  the 

foliage,  as  in  an  average  or  effective  medium.  This  medium  is  characterized  by  the  relative 

dielectric  constant  cr  where 
P 

trp  =  kplko  p  =  vorh  (B\) 

When  these  coherent  fields  strike  the  foliage-air  interface,  there  will  be  reflected  and  transmitted 
fields.  The  reflected  field  is  ignored  because  the  dielectric  contrast  between  free  space  and  the 
average  foliated  medium  is  not  very  significant.  However,  a  proper  accounting  for  the  change 
in  medium  must  be  undertaken  for  the  transmitted  field. 

The  average  incident  field  at  z=h,  coming  up  through  the  foliage,  is  given  by  (1.88),  i.e., 

<  Zsfl  -  Sw  +  LFKzf, L  m 

The  k  superscript  on  the  fields  and  the  k  subscript  on  the  operator  K £  indicate  that  the  medium 
is  characterized  by  the  wavenumber  associated  with  the  average  field  in  the  foliage  medium. 
We  now  have  to  continue  this  field  into  free  space  above  the  foliage.  It  is  well  known*  that  the 
fields  above  the  z  =  h  plane  can  be  determined  from  the  spectral  representation 

<  E  (rt,  2)]  =  (2 n)~2f  F{kt)  exp  (  -;'k[z  -  /i]  -jkt .  r()dkt  (M) 

where 

{kl-krkt)XI2 

K  = 

-j(kt .  kt  - 

and  dkt  =  clkxdk„'  kt  =  k^x  +  kyy .  liquation  (B3)  represents  the  average  field  above  the 
foliage  as  a  superposition  of  plane  waves  whose  complex  vector  amplitudes  arc  given  by  l'[kt). 

*  Collin.  R.F.,  Antennas  and  Radiowave  Propagation ,  McGraw-Hill,  New  York,  1  ‘ ? 8 5 . 


k0>kr  kt 


kt  •  kt  >  k~ 


m 
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Since  <  E  (rt  z)]  satisfies  Maxwell's  source  free  equations,  we  must  have  7  .  <  E  (rt,  z)]  =  0. 
This  implies  that  all  the  vector  components  of  F  arc  not  independent  and,  in  fact, 

(kt  +  kz)  •  F  =  0 

or 

(V*  +  kyFv) 

fz - .  —  m 

The  x  and  y-components  of  F ,  denoted  as 

Ft-Fjc  +  /y,  (56) 

may  be  determined  from  (B3)  by  setting  z=h  and  taking  the  inverse  two-dimensional  Fourier 
transform  of  both  sides  to  yield 

Ft(k{)  =  /  <  Et(rt,  z  =  A)]  exp  jk( .  rt)drt  (Bl) 

In  (B7) 

<  £f(r,,z  =  A)]  =  <  Ex(rrz-h)]x  +  <  Ey(rt ,  z  =  A)]y  (58) 

is  taken  to  be  the  x  and  y  components  of  the  field  in  (B2);  that  is,  the  upward  traveling  field  in 
the  foliage  layer.  Thus,  (B7)  can  be  written  as  follows; 

F((kt)  =  /  <  Esj(r(<  z  =  A)]  exp  jkt .  r[)dr[  (59) 

When  the  point  of  observation  is  sufficiently  far  from  the  illuminated  area  on  the  interface 
at  z-  h,  the  integral  in  (B3)  can  be  asymptotically  evaluated  to  yield 

<  £  (rf,  z  >  >  A)]  ~ j2k0  cos  9SG(R)F(  -  kS{)  (5 10) 

where 

(7(5)  =  exp  (  —  jkQR)IAnR  (511) 

and 

kS(  =  —  kQ  sin  0S  cos  —  k()  sin  0S  sin  4>^v  (512) 
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The  angles  8S  and  4>s  are  the  polar  and  azimuthal  scattering  angles  while  R  is  the  distance  from 
the  top  of  the  foliage  layer  to  the  point  of  observation.  In  this  limit,  all  we  need  is  the 
two-dimensional  Fourier  transform  of  the  total  average  field  at  z  =  h.  In  terms  of  the  parts  o( 
<  Esj\,  as  given  by  (B2),  the  overall  process  of  getting  from  the  current  on  the  rough  surface  up 
through  the  foliage  and  into  free  space  is  a  bit  more  involved. 
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Appendix  C 

The  Conversion  of  Fields  Into  Scattering  Cross  Sections 

While  the  fields  obtained  in  the  main  body  of  this  section  are  fundamental  quantities  they 
are  not  the  most  useful  to  a  radar  engineer.  For  incoherent  radars  and  even  most  coherent 
radars,  the  characteristic  of  primary  interest  is  the  scattering  cross  section  of  the  terrain.  This 
quantity  is  related  to  the  second  moment  of  the  scattered  field  so  it  is  determined  by  the 
fluctuating  component  of  the  scattered  field  in  all  directions  except  the  specular  direction.  In 
the  specular  direction,  the  scattered  power  comprises  both  coherent  and  incoherent  parts 
although  the  coherent  part  is  frequently  negligable.  If  the  coherent  power  is  not  ncgligable,  it 
gives  rise  to  a  range-dependent  scattering  cross  section. 

With  scattering  from  bare  ground,  the  cross  section  of  interest  is  the  scattering  cross 
section  of  the  surface  per  unit  scattering  area  or  a°  .  This  normalized  monostatic  cross  section 
is  determined  from  the  fluctuating  part  of  the  scattered  field  as  follows; 

a°pq  =  Um  {4 nR2  <  \  6ES .  q  \ 2  >  / 1  £,  | 2 A)  (Cl) 

A—*  oo 

where  R  is  the  distance  from  the  radar  to  the  midpoint  of  the  illuminated  area,  A  is  the 
illuminated  area  on  the  surface,  and  £,•  is  the  incident  field.  The  subscripts  p  and  q  denote  that 
the  polarization  p  is  transmitted  while  the  q-component  of  the  backscattered  field  is  sampled. 
This  normalized  cross  section  depends  on  only  two  radar  parameters,  namely,  the  frequency  and 
the  polarization. 

The  normalized  cross  section  a °  is  adequate  for  bare  surfaces  and  surfaces  having 
vegetation  layers  which  are  small  in  depth  compared  to  the  range  resolution  of  the  radar. 
However,  when  the  foliage  depth  exceeds  the  radar  range  resolution,  the  surface  scattering  cross 
section  should  be  replaced  by  a  volume  scattering  cross  section.  That  is,  instead  of  (Cl),  the 
following  volume  cross  section  should  be  used; 

<7*  =lim  {4nR2<\dEs.q\2>l\Ei\2l'}  (C2) 

V— *oo 

where  V  is  the  illuminated  volume  and  R  is  the  distance  from  the  radar  to  the  midpoint  of  the 
volume.  It  should  be  noted  that  just  as  (Cl)  assumes  completely  incoherent  surface  scatter  so 
docs  (C2)  assume  completely  incoherent  volume  scatter.  That  is,  under  the  complete 
incoherence  condition,  cither  a 0  pqA  or  apq V  (whichever  is  appropriate)  can  be  used  in  the 

conventional  radar  equation  for  the  target  scattering  cross  section  to  give  the  power 
backscattered  by  the  terrain. 
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As  summarized  in  Section  1.4.1  there  are  essentially  three  fields  to  consider  when  dealing 
with  the  scattering  from  vegetation  covered  terrain.  The  first  field  is  due  to  scattering  by  the 
foliage  of  a  field  incident  from  free  space.  This  field  is  a  fluctuating  field  in  that  the  foliage  layer 
is  assumed  to  be  so  large  in  terms  of  a  wavelength  that  there  is  essentially  no  coherent 
backscattering.  The  other  two  fields  are  the  coherent  and  incoherent  fields  w'hich  have 
propagated  down  through  the  foliage  and  have  been  scattered  by  the  surface  back  up  through 
the  foliage.  Of  these  three  fields,  the  two  of  most  interest  are  the  fluctuating  field  scattered  by 
the  foliage  on  its  way  to  the  surface  and  the  fluctuating  field  scattered  by  the  surface  back  up 
through  the  foliage. 

Scattering  of  the  free-space  incidence  field  (propagating  down  to  the  surface)  by  the  foliage 
is  described  within  the  distorted  Born  approximation  by  (1.86).  The  mean  square  value  of  this 
quantity  may  be  written  as  follows; 

<\dEf(j\‘'>  =  <  LK-£  <  Ej-j]  .  LAV  <  £yy]  >  -  <  LA^  <  Ej-j]  >  <  LAV  <  Zyy]  >  (C3) 

If  the  point  of  observation  is  sufficiently  far  removed  from  the  scattering  region,  the  far  field 
approximation  for  the  Green's  function  in  (1.11)  may  be  used  to  yield 

klsiR)  £  (er„  ~  [)sn  exp[/*<A  •  (^  +  ^)]{  1  -  ks(ks  • }  (C4) 

n=\ 

A 

where  ks  denotes  the  direction  of  the  scattering  observation.  Rewriting  (C4)  as 

=  2  N  -  - 
KZ  =  koi(R)  £  Qn(ra  +  rn) 
n=  1 


permits  the  following  expression  for  the  first  term  on  the  right  side  of  (C3); 


LKZ  <  EfS  •  LKZ  <  Rfd\  =  (7^-) 2 ^  I LQn(T0  +7n)<  ~EfJr0  +  ^)] 1 2 


Ko  ,2 


N  _  _ 


N 


+  (  4_  d  )  —  + 'n)  <  (r0  +  rfj)]  •  —  kQm{r0  +  rm)  <  Efd^o  +  rm^  (C5) 

H7rn  rt=l  m- 1 


m^n 


What  (C 5)  docs  is  to  separate  out  the  N  product  terms  which  contain  like  elements.  Assuming 
that  the  scattcrcrs  arc  independent  of  each  other,  an  assumption  which  is  consistent  with  the 
distorted  Born  approximation,  an  average  of(C5)  yields 


<  LAV  <  Lyj]  •  LAV  <  F.fj] 


4  nR 


■f  -V  /  <  I  LQn(r0  +7n)<  Ejd  (rQ  +  rn)]  | 
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+  (-^-)2'v2(!  -  llN)\fL  <  Qn(r0  +7*)]  <  Ef(f0  +7n)]p(rn)drn\2  (C6) 

where  p(rn)  is  the  probability  density  function  for  the  location  of  the  scattcrers,  symbolized  by 
the  vectors  rn  n-  1,2, As  N  -*  oo,  the  second  term  on  the  right  side  of  (C6)  goes  to  the 
negative  of  the  second  term  on  the  right  side  of  (C3)  so  they  cancel  and  this  leaves 

k2  _  _ 

<  1 8Efd\ 2  >  =  J  <  !  LQn(k0  +7n)<Ejd(70+  7ny]  1 2Mr„)d7n  {Cl) 

The  < .]  brackets  inside  the  integrand  denote  a  conditional  average  over  scatterer  size,  shape, 
orientation,  and  dielectic  constant. 

Substituting  from  (C4)  and  simplifying,  yields  the  following  result; 

<\6Efd\2>  =  {-~-)2  N  <  (irn  -  1  )2  >  //  /  [  1  -  ( V  e/d)2]  <  Efd (70  +7n)]<  E}d %  +  7n)] 

<VS> 

.  p(7n)  exp[ fkjts  •  (r0  -  7'0  )]dr0dr  'Q  drn  (C8) 


where 


<  £/</]  =  <  Efcti  Jd 

and  <  Vs>  is  the  average  volume  of  a  scatterer  and  shape  of  a  scatterer.  For  plane  wave 
illumination  of  the  foliage  layer,  the  average  field  in  the  layer  is  given  by 

<  E[(f0  +  Ci>]  =  Ei  exPC  ~iki(  •  (fot  +  rnt)  ~ )kp70  +  ’«)]  (C9) 


where 


k:  =  k:  X 
H  ‘x 


t  A 

+ v 


(CIO) 


with  kj  and  k(  given  by  (1.5 la)  and  (1.51b),  and  kp  is  the  propagation  constant  of  the  average 
medium  given  by  (1.54).  The  average  field  has  its  transverse  (to  z)  wavenumber  equal  to  k(  in 
order  to  be  in  phase  synchronism  with  the  free  space  incidence  field  across  the  upper  average 
boundary  (z  =  h).  Substituting  (C9)  in  (C8)  yields 


1 8Ej-d\ "  >  =  (■ 


4  nR 


•)2  .V  <  (irn  -  1 7>  ff  f  [!-(*,-  e/dT]  I  C/I  *  exp  \J{kS{  -  .  (r„(  -  r  ) 

<»/5> 
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where  the  conjugate  on  kp  =  kg  -  jk‘p  is  necessary  because  the  propogation  constant  for  the 
average  field  may  be  complex.  (Cl  1)  may  be  rewritten  as  follows, 

,2. 

<  \6Efd\2>  =  (-^)2  N  <  (crn-\)2>  |£,|2/l  /  [l-(4*^>2] 

<vs> 

•  exp [j(kSt  -  kQ  .70(  +j(kSz  -  kp)z0]d70 1 2  exp  (  -2kgz„)p(rn)drn  (Cl 2) 

The  integral  over  <  Vs  >  has  the  same  form  as  the  scattering  pattern  of  a  scatterer  of  volume 
<  Vs>  and  supporting  an  internal  field  of  the  form  in  (C9).  Thus,  with 

P(kjt  +  kpz;  ks)=  f  U-(ks‘ej-dfye\p[J(kst- k^.70t+j(kSz- 'ip)Zo]d70  (C13) 

<y,> 

(Cl 2)  becomes 

k2  _ 

< \6Efd\2>  =  (-^-)2N<(irn-D2)>  \Ei\2\P(ki(  +  k^,ks)\2  (C14) 


.  /  exp  (  -2 klp?n)p(rn)<frn 


If  the  scatterers  are  assumed  to  be  uniformly  distributed  about  the  volume  V  -  A  .A  where  h  is 
the  depth  of  the  foliage  layer  and  A  is  the  area  in  the  z  =  0  plane  (sec  Figure  1.3)  then 


and  (04)  becomes 


PCn)  = 


1  IAh 


Jn  in  V 
rn  not  in  V 


(Cl  5) 


ko  ,2 


<  \6Efd\Z>  =  (^PV  <(crn-\r)>  |£,M  P(kit  +  kpz ;  ks)  l 


[1  -cxp(_2  k'ph)] 


Substituting  this  into  (C2)  yields  the  following  result  for  the  scattering  cross  section  per  unit 
scattering  volume; 


(Cl  6) 


°V  =  ~v~£n  <  (' %  ~  1)2  >  I  ^  +  kf ;  A^)  1 2[1  -  exp(  -  2A£/j)]/2A^A 
If  the  scattering  amplitude  is  normalized  to  its  maximum  value,  i.e.  <  Vs>  ,  (C16)  becomes 

*■ 4  ,  -  -  , 

"V~<  ^><(tr/|-l)S|P(/t(-+yA;^)/<  r 

.vc[l-exp(-2^)]/(2^A)  (C17) 

where  vc  is  the  average  fractional  volume  occupied  by  the  scatterers  or 

vc  =  N<Vs>IV  (C18) 

This  result  differs  from  what  would  be  obtained  with  previous  results  in  the  appearance  of  the 
factor  2k‘ph.  That  is,  previous  results  would  erroneously  replace  this  by  4 kj/i. 

A  consequence  of  the  sparse  scatterer  assumption  is  that  (Cl 7)  very  nearly  splits  into  the 
product  of  two  distinctly  different  factors.  The  first  factor,  with  the  exception  of  its  dependence 
on  the  propagation  constant  of  the  average  medium,  is  a  function  of  the  properties  of  a  single 
average  scatterer  only.  The  other  factor  depends  on  average  medium  parameters  only.  Since  P 
or  the  scattering  pattern  is  the  only  function  in  (Cl 7)  which  depends  on  the  scattering  direction, 
this  shows  that  the  "scattering  pattern"  of  the  ensemble  is  determined  largely  by  the  scattering 
pattern  of  an  average  scatterer. 

The  computation  of  the  normalized  scattering  cross  section  due  to  the  fields  which 
propagate  down  through  the  foliage  strike  the  underlying  surface,  and  then  are  rescalicred  back 
up  through  the  foliage  is  considerably  more  difficult.  This  is  due  primarily  to  the  fact  that  there 
are  essentially  three  random  scattering  processes  occuring  in  this  situation.  Equation  ( 1 .87) 
gives  the  weak  volume  scattering  approximation  for  the  fluctuating  field  described  above,  e.g. 

4/  -  4u  +  ^(1-  +  LpKj$su  +  4«>  (C19) 

Eftsu  is  the  field  scattered  into  free  space  by  the  rough  surface  when  illuminated  by  SEjj.  The 
term  L(  1  -  P)K^E^SU  represents  the  fluctuating  field  scattered  by  the  particles  when  supporting 
the  field  £(i  E^u  is  the  field  scattered  into  the  average  medium  by  the  surface  when 
illuminated  by  <  Ejj>.  The  term  L(  \  -  /’)AV/:£5U  is  the  fluctuating  field  scattered  by  the  particles 
when  supporting  the  field  E^u  .  Finally,  the  term 

L(1  -  r)KzLPKZkE^u 
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represents  a  double  scattering  of  the  fluctuating  field  E$su  and  it  is  not  significant  compared  to 
the  other  terms  in  (C19)  when  there  is  a  sparse  population  of  scatterers.  Thus,  (Cl 9)  can  be 
reduced  to  the  following; 


6Esf  *  ESsu  +  L(  1  —  P)Ez{EKasu  +  ESsu)  (C20) 

Assuming  that  there  is  no  correlation  between  these  various  fields  and  recognizing  that 
each  are  zero  mean  quantities  results  in  the  following  expression  for  the  mean  square  value  of 

SEsf  '• 

<  \SEs/\2>  =  <\E5su\2>  +  <  \L{\-P)KIEkasu\2> 

+  <  \L(l-P)KrE6su\2>  (C21) 

The  first  term  on  the  right  of  (C21)  is  due  to  surface  scattering.  After  considerable  algebra,  it 
may  be  reduced  to  the  following  form; 

<  i 2 > 

•  exp  [jk0ks  •  (rs  .7;  )]  >  dr^  dr ^  (C22) 

where  *\0  is  the  impedance  of  free  space,  is  the  current  induced  on  the  rough  surface  by  the 
downward  propagating  foliage  field  <5£y^  (and  satisfies  equation  (1.64), 

dns  =  dxsdys  dr'ts  =  dx*  dy'z  , 

A  — * 

and  denotes  the  direction  of  J ^  This  result  can  be  reduced  to  a  surface  scattering  cross 
section  through  the  use  of  (Cl),  e.g. 

k2  2 

-  (— <  [1  -  &s-hj)2V5fs)Jdfis  ) 

4n\Ej\  / 


.  e\p[jk()ks .  (rs  .  r  '  )]  >  drtdr  ^  (C22) 

This  is  about  as  far  as  this  formulation  can  be  taken  without  a  specific  form  for  the  current. 
However,  it  must  be  remembered  that  the  average  in  (C23)  involves  an  average  over  the  random 
properties  of  the  surface  and  the  random  properties  of  the  discrete  scatterers  including  their 
positions.  The  .l-*  factor  in  (C23)  will  be  cancelled  by  one  of  the  integrations  over  y  and  y. 
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For  the  second  and  third  terms  in  (C20),  the  following  relationship  developed  in  an  earlier 
part  of  this  appendix  is  useful; 


< \L{\-P)KLE\2>  =(-^-)2<(^-i  )2>vf 


/  <[!-(*, 


A .  2-* 

e)  ] 


Vs> 


•  E(ro  +  rn>E  Oo  +Zi)>  exP[/Vi  *  (r0  - ~r'o  )1drod^o  (C23) 

The  averaging  in  (C23)  inside  the  integral  is  over  the  random  positioning  of  the  scatterers  (for 
the  upward  traveling  field),  the  random  surface,  and  the  random  locations  of  the  scatterers  for 
the  downward  traveling  field.  The  form  of  (C23)  suggests  that  to  a  good  first  approximation 
the  second  and  third  terms  in  (C20)  give  use  to  a  scattering  cross  section  per  unit  scattering 
volume  similar  to  the  one  in  (Cl 6)  but  modified  by  the  random  surface  effects.  For  example, 
the  second  term  on  the  right  of  (C20)  gives 

k 4 

°Vasux  47  <  Vs>K  (ern  ~  l)2  >  I  P(ks,  + kf-  ks)l  <  Vs>  1 2 


.  (vc[l  -  exp(  - 2&^A)] j(2k‘ph)} Sasu 


(C24) 


while  the  third  term  yields 


V 

°6su  ~ 


vs><brn-\)2>\p£st  +  k*\ 


ks)l<Vs>\2 


exp(  -2klD 


osu 


where 


^•/i  C)  /*  i*  A  A  9 

Sasu  =  —  If  <  [1  -  (*,  ’Ja/VajCs)J\jir;  ) 

•  exp [jk0kst  •  (r*tf  )  +jkp(C  -  C')]  >  dTtdT ^ 

^  2  2 

''i>SU  =  If  <  [1  -  (ks»jty)  ) 


•  exp[ /Vj  •  (fj  -  T;  )]  >  drsdrlt 


(C25) 


(C26) 


(C27) 
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Of  course,  the  averaging  in  (C26)  and  (C27)  must  include  not  only  the  random  surface  but  the 
discrete  scatterers  as  the  field  is  scattered  down  to  the  surface. 

To  the  level  that  has  been  developed  above,  the  effects  of  a  random  foliage  layer  on  a 
random  surface  may  be  described  by  a  scattering  cross  section  and  two  volume  cross  sections. 
It  should  be  remembered  that  no  coherent  fields  reflected  from  the  surface  have  been  included. 
Clearly,  there  is  much  work  to  be  done  on  this  problem.  For  example,  one  of  the  things  that 
should  be  done  next  is  to  use  some  classical  approximations  for  the  surface  current  (such  as 
physical  optics  and  perturbation)  and  carry  the  above  calculations  through  to  completion.  This 
would  certainly  give  a  better  feel  for  the  relative  importance  of  the  various  effects  above. 


2.0  Asymptotic  Surface  Scattering 

In  Section  1.3  we  reviewed  a  number  of  approximate  methods  for  dealing  with  surface 
scattering  and  their  relationship  to  iterations  of  the  MFIE.  Using  these  iterates  along  with 
asymptotic  techniques  to  determine  the  scattered  field  provided  further  insight  into  the 
limitations  and  capabilities  of  the  iterative  approach.  However,  the  iterative  approach  is  not 
necessarily  the  one  to  use  to  improve  on  known  asymptotic  results  because  such  improvements 
may  possibly  comprise  the  partial  summation  of  an  infinite  number  of  iterates.  There  are  other 
methods  for  extending  asymptotic  results  and  the  purpose  of  this  section  is  to  apply  and  develop 
one  such  technique.  These  asymptotic  extention  methods  generally  augment  iterative  solutions. 

2.1  Luneburg  -  Kline  Expansion  for  the  Surface  Current 

2.1.1  Introduction 

One  of  the  few  scattering  properties  that  are  known  about  certain  classes  of  rough  surfaces 
is  how  they  scatter  electromagnetic  energy  in  the  high  and  low  frequency  limits.  The  caveat 
"certain  classes"  is  necessary  because  even  in  these  asymptotic  frequency  ranges,  there  arc  some 
surfaces  which  are  not  amenable  to  classical  analytical  techniques.  For  these  surfaces  or  more 
general  ones,  there  are  problems  with  existing  theories  or  models  when  the  frequency  of  interest 
is  between  the  low  or  high  asymptotic  limits  (!].  The  problems  are  not  with  a  lack  of  models 
but  with  a  good  understanding  of  the  limitations  of  the  models.  For  example,  there  has  recently 
been  some  very  good  work  appear  on  a  phase  perturbation  approach  to  rough  surface  scattering 
[2,  3].  This  Rytov-like  method  has  shown  promise  of  extending  the  range  of  conventional  power 
series-like  Held  perturbation  theory  to  the  point  where  the  rms  surface  roughness  is  the  order 
of  an  electromagnetic  wavelength  [4].  Of  equal  importance,  however,  is  the  fact  that  some  of 
the  advocates  of  this  approach  have  very  carefully  studied  the  validity  of  this  method  by 
applying  it  to  deterministic  surfaces.  This  type  of  validation  analysis  must  be  done  if  the 
method  is  to  be  useful  as  an  engineering  tool. 

Phase  perturbation  is  essentially  a  low  frequency  method  whose  upper  frequency  range  has 
been  extended  by  partial  summation  [3j.  This  approach  of  extending  a  low  frequency  method 
into  higher  frequency  ranges  is  fairly  common.  It  is  based  on  the  fact  that  most  low  frequency 
techniques  do  have  higher  order  "terms"  which  can  formally  be  derived;  some  of  these  terms  can 
also  be  computed.  The  problem  with  such  extensions  is  that  it  is  never  quite  clear  when  the 
higher  terms"  converge  and  when  they  give  meaningless  results.  Little,  if  any.  work  has  been 
done  in  the  rough  surface  literature  on  extending  high  frequency  solutions  into  lower  frequence 
regions.  However,  a  basic  technique  for  doing  so  has  existed  for  a  number  of  years;  that  is.  the 
Luneburg- Kline  (LK)  representation  for  the  field  scattered  by  a  large  object  [5j.  The  classical 
Luneburg-Klinc  (LK)  representation  is  a  scries  in  powers  of  the  electromagnetic  wavelength  (or 
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inverse  powers  of  the  electromagnetic  wavenumber)  for  the  field  scattered  from  a  body.  The 
zeroth  order  or  wavelength  independent  term  corresponds  to  the  geometrical  or  ray  optics 
prediction  for  the  reflected  field.  This  term  depends  on  the  curvature  and  the  reflection 
coefficient  of  the  surface  at  the  point  of  reflection  (the  so-called  specular  point).  The  term  that 
varies  as  k~ *  also  depends  on  these  same  quantities  but,  in  addition,  appears  to  oe  sensitive  to 
the  distance  from  the  specular  point  to  the  shadow  boundary  [6].  This  distance  dependence 
indicates  the  nonlocal  nature  of  the  first  wavelength  dependent  term. 

Lee  [7]  has  recently  obtained  the  k~ 1  coefficient  for  the  scattered  field  and  the  current  on 
a  perfectly  conducting  convex  body.  This  was  a  significant  result  but  it  was  a  bit  too  restrictive 
for  application  to  the  rough  surface  problem.  For  example,  Lee's  solution  must  be  augmented 
by  possible  multiple  reflections  on  the  surface,  the  presence  of  a  creeping  wave  where  the 
incident  field  just  grazes  the  surface,  and  any  possible  edge  diffraction  where  the  surface  has  a 
relatively  sharp  edge-like  behavior.  In  addition,  the  solution  obtained  by  Lee  was  very' 
complicated  from  an  algebraic  point  of  view  and  involved  surface  characteristics  whose  statistics 
are  not  necessarily  w'ell  know. 

Thus,  with  an  eye  toward  possibly  generalizing  Lee  s  results  and  also  obtaining  more 
insight  into  the  L-K  approach,  an  alternate  methodology  was  developed  IS].  First,  rather  than 
expanding  the  scattered  and  incident  fields  in  L-K  series  and  then  applying  the  boundary 
conditions  to  obtain  the  expansion  coefficients,  the  surface  current  induced  on  the  surface  by 
the  incident  field  was  expanded  in  an  L-K  series.  This  expansion  was  used  in  the  Magnetic  Field 
Integral  Equation  (MFIE)  to  generate  a  hierarchy  of  integral  equations  for  the  expansion 
coefficients.  This  sequence  of  integral  equations  exhibited  some  rather  interesting  properties. 
First,  it  was  recurrsive  in  that  knowledge  of  the  nth  integral  equation  solution  determined  the 
n+  1  integral  equation  solution.  Second,  the  "source"'  or  Born  term  in  these  integral  equations 
depended  entirely  on  the  asymptotic  evaluation  of  a  known  integral  in  inverse  powers  of  k0. 
Finally,  and  most  interesting  of  all  was  the  fact  that  each  integral  equation  could  be  solved  as 
they  were  all  just  like  the  equation  obtained  by  iteratively  solving  the  MFIE  in  the  high 
frequency  limit. 

There  were  two  immediate  consequences  of  this  work.  First,  it  was  found  that  the  L-K 
representation  for  the  current,  even  though  forced  to  satisfy  an  exact  integral  equation, 
produced  a  zero  current  on  the  part  of  the  surface  shadowed  from  the  incident  field.  That  is. 
even  the  higher  order  terms  in  \,  k0  were  identically  zero  in  the  shadow  zones  of  the  surface. 
Although  it  is  not  directly  obvious  why  this  is  the  case,  it  is  apparent  that  since  the  L-K 
representation  is  asymptotic,  it  fails  to  converge  to  the  value  of  current  which  is  correct  for  ail 
frequencies.  Thus,  in  the  shadow  zones  of  the  surface,  the  L-K  representation  is  valid  in  the 
optical  limit  only.  This  result,  in  itself,  is  significant  because  it  is  the  first  time  that  a  limitation 
of  the  L-K  representation  has  been  obtained  via  an  exact  analysis.  A  second  major  result  of  this 
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work  was  to  show  that  all  the  L-K  expansion  coefficients  could  be  determined  from  the  behavior 
(as  a  function  of  k0  )  of  a  two-dimensional  integral  of  known  functions.  Although  time  did  not 
permit  the  evaluation  of  this  integral,  there  are  techniques  for  doing  so  [9]  and  it  is  suggested 
that  this  is  a  useful  avenue  for  future  work. 

2.1.2.  Analysis 

The  problem  to  be  addressed  here  is  the  determination  of  the  electric  current  density  Js 
induced  on  the  surface  z  =  C(rJ)  by  an  incident  magnetic  field  H‘.  The  surface  z  =  £(rf)  separates 

A 

free  space  (z>C)  from  a  perfectly  conducting  medium  (z<,  C).  The  unit  vector  k-t  specifies  the 
direction  of  travel  of  the  incident  field.  The  electric  surface  current  density  Js  must  satisfy  the 
Magnetic  Field  Integral  Equation  (MFIE)  as  follows; 

Js(r)  =  2n(r)xHl(r)  +  2n(r)xfJs(r0)xV0G(r  -  r0)dsQ  (2.1) 

In  (2.1),  n(r )  is  the  unit  normal  to  the  surface  at  the  point  r  =  rt  +  C(rf)z  and  is  given  by 

n(7)  =  [  -  t;J  -  Cyy  +  z]/(l  +  cl  +  Cp 1/2  (2.2) 

where  =  d£ldx  ?nd  t 'y  =  d^dy  are  the  x  and  y  surface  slopes  at  the  point  r  on  the  surface. 
G(r  -  r0)  is  the  free  space  Green's  function,  i.e. 

G(r  -  r0)  =  exp(  -jkQ\7  -  rQ |  )/4rr  |7  -  7a\ ,  (2.3) 

and  V0  is  the  conventional  three-dimensional  gradient  evaluated  on  the  surface  (z  =  £  and 
zo  ~  C o  )•  Noting  that  the  area  integration  over  the  surface  can  be  converted  to  one  over  the 
z  =  0  plane  through 


ds0  =  (i  +  c2  +  $1/2  J7lo 

where  dr(°  =  dx0dy0,  (2.1)  can  be  rewritten  as  follows; 

J(r)  =  2.V  xll‘  +  2.V  xfj  {rQ)xl 0G{  \r  -  r0\)dr,^ 


(2.4) 


where 


•/('•)  =  J%(r  )(1 


(2.5) 


and 
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A 

2 


(2.6) 


/V  =  -tj  -  tyy  + 

Ify  (r )  can  be  determined  from  (2.4),  the  scattered  field  can  be  found  from  the  following  integral 
expression; 


HS(R )  =  7  x/7(T)G(  |  R  -  7 1  )dTt  (2.7) 

As  a  preparatory  step  to  introducing  the  Luneburg- Kline  expansion,  the  k0  dependence 
introduced  by  the  incident  field  is  removed.  That  is,  with 

J  {r  )  =  L  (r  )  exp  ( -jk-t  .  r  ), 

H\r)  =  ^  exp  (~jkr7), 

and  ^  a  constant,  (2.4)  may  be  rewritten  as  follows; 

L(r)  =  2N  x  ^  +  2/V  x  /  L(r0)x?0G(|r  -r0l)  exp  (jki  •  Ar  )drt°  (2-10) 

where  Ar  =  r  -  The  purpose  of  (2.8)  is  to  remove  the  known  high  frequency  behavior  or 
dependence  on  k0  from  the  current.  That  is,  as  k0-> oo  it  is  known  that  L(r)  is  independent  of 
kQ.  The  modified  current  L  (r  )  is  expanded  in  a  Luneburg- Kline  series,  i.c. 


(2.8) 

(2.9) 


L(r) 


v  jf'^) 
n=o  k"o  ’ 


(2.11) 


where  the  vector  expansion  coefficients  jn(r  )  ,  n 38  0,1,...,  are  independent  of  the  electromagnetic 
wavenumber  k0.  (2.11)  is  next  substituted  into  (2.10)  and  it  is  assumed  that  term  by  term 
integration  is  permissible*  so  that  the  following  results; 


X  Ur  )k0  n  =  2<V  xhj  +  2,V  x  T  kQ  n  f  jn{r0)  xV0G  exp  (jkj  •  Ar  )drt  (2.12) 

n= 0  n= 0  0 

In  order  for  (2.12)  to  be  satisfied,  the  integral  term  must  also  have  a  Luncburg-Kline  expansion. 
That  is,  if 


*  This  assumption  is  tantamount  to  requiring  that  the  L-K  scries  be  uniformly  convergent  over 
r, ^  (  -oo.co);  a  requirement  that  is  probably  not  satisfied. 
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(2.13) 


ft(ft  ft)  =  f  jn(r0)  xl  QG  exp  (ft  •  AT)cfrft 
then  it  may  be  written  as  follows; 

ft(ftft)  =  X  Tnm(?)lk™  (2.14) 

m=0 

where  the  7ft,,  expansion  coefficients  are  independent  of  k0.  Combining  (2.13)  and  (2.14)  yields 

Tno  +  'TL  +  - T-  +  -  =  /  x  7  0G(  |  aT  | )  exp  (ft  •  \r)dr,  (2- 1 5) 

k°  ft 

so  that  the  vector  expansion  coefficients  can  be  determined  as  follows; 

ft0  =  lim  f  jn(r0)  x  Vft7(  |a7|)  exp  (ft  .  Ar  )<ft  .  (2. 16a) 

ft—co  w 

7, 1  =  Um  ft  [  /  jn  X  7  0G  exp  (ft  •  ^r)drt^  -  Tno  1  (2. 1 66) 

AC^  *00  L 


rwi  ft 111X1  ft 
ft-*00 


-*•  —  _»  m~  1  ftp 

f  JnxloG  exP  (ft  • Ar  )drta  -  \ 

p=0  ft 


(2.16c) 


Substituting  (2.15)  into  (2.12)  yields 


I  ft'  )ft  =  2.V  x 

«=0 


6/  +  S  S  ftm(r)ft 

rt=0  A72=0 


—n—m 


(2.17) 


so  that  equating  like  powers  of  ft  yields 


jjr)  =  2-v  x  ft-  +  ft0) 

(2.18a) 

7|(r  )  =  2.\  x  f  /  |(j(r  )  +  T ()|(r  )} 

(2.1 86) 

y2(r  )  -  2.V  x(Tn(r)  +  'ftylr  )  +  7'()2(c  )} 

(2.18c) 

etc. 
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Substituting  from  (2.16)  into  (2.18)  yields  the  following  sequence  of  integral  equations  for  the 
vector  expansion  coefficients; 

jJr)  =  2  N  x  ht  +  2jV  x  lim  f  j0(r0)  *  V0G(\Ar  I )  exp  (Jkj .  Ar  )dr,  (2.19a) 

Kq—*oo  0 

j\(r)  =  2 N  x  lim  (k0  [  f  j0(r0)  x  VflG(fAr  I )  exp  (jkj  .  A r  )drt  -  TQ  J} 

kg—^OO  L  O  J 

4*  2 iV  x  lim  / y’i(r  )  x  V0G(|  Ar  | )  exp  (jk:  •  Ar  )drt  (2.196) 

^o°  o 

•  • 

•  • 

•  • 


This  sequence  of  integral  equations  has  a  number  of  interesting  properties.  First,  except 
for  the  source  or  Born  term,  all  the  integral  equations  are  identical  in  form  in  that  they  appear 
as 


jp{r  )  =  2/V  x  7p(r  )  +  2/V  x  lim  f  jp(r0)  x  1 0G{  I  Ar  I )  exp  (jkj .  Ar)drt^ 


(2.20) 


p  =  0,1,2,... 

Also,  since  the  source  term,  2N  x  ~sp{r),  in  (2.20)  depends  on  j0,  y'j,...,  and  jp_ j,  the  integral 
equations  are  recurrsive.  Thus,  if  j0  can  be  determined  then  it  should  be  possible  to  determine 
all  higher  order  vector  coefficients.  Because  the  vector  expansion  coefficients,  jn(r ),  are 
independent  of  k0,  it  is  further  noted  that  the  source  term  sp(r  )  is  determined  almost  completely 
by  a  Luneburg-Kline  series  representation  for  the  following  integral 

/  V,/7(|a7|)  exp  (Jkj  •  Ar  )drt  . 

For  example,  from  (2.19a)  it  follows  that 


s()(r  )  =  hj 


(2.21a) 
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while  from  (2.19b) 


7[  (r  )  =  lim  k0\  f  j0(rQ)  x  7  flG(  |  a7  | )  exp  (/*• .  a7  )drt 

kQ — *00  L  ° 

(2.216) 

-  ,  lim  /  /0(70)  x  7  G(  I  aT  | )  exp  (/*f-  •  a7)c/7> 
k0^oo  o  j 


It  should  be  noted  that  the  presence  of  y0(r0)  under  the  integral  signs  in  (2.21b)  has  a  very 
minimal  effect  because  the  dominant  terms  are  the  ones  which  depend  upon  k0. 

One  final  but  very  important  point  about  the  sequence  of  integral  equations  in  (2.19)  is 
that  their  solutions  are  known.  In  fact,  (2.19a)  is  a  slightly  altered  form  of  the  Magnetic  Field 
Integral  Equation  in  the  high  frequency  (kQ^oo)  limit,  so  its  solution  is  given  by  [10] 

2j V  x  hj  (r  not  shadowed) 

Jo?)  ~ 

0  ( r  in  shadow) 


There  is  also  the  possibility  of  a  multiple  scattering  contribution  to  ja{r )  from  other  points  on 
the  surface  [10].  Thus,  for  the  nth  vector  expansion  coefficient,  the  solution  is  (except  for  the 
contributions  of  multiple  ray  bounces  on  the  surface) 

2 N  x  sn{r  )  (r  not  shadowed) 

7nt )  -  (2-22) 

0  (r  in  shadow) 


The  complete  surface  current  density  is  thus  obtained  by  combining  (2.5),  (2.8).  (2.11),  and 
(2.22),  i.e. 


7jr)  -  0  +  <£  +  C h  ll2exp(-jkr7)  I  (2.23) 

n= 0  kQ 

One  of  the  immediate  consequences  of  (2.23)  is  that  the  current  on  any  shadowed  portion 
of  the  surface  is  identically  zero.  This  is  obviously  a  high  frequency  approximation,  but  the 
analysis  presented  above  makes  no  explicit  approximations  and,  in  fact,  appears  to  be  exact. 
Clearly,  an  exact  analysis  cannot  lead  to  an  approximate  result.  What  is  happening  in  this  case 
is  that  the  L-K  scries  is  doing  the  best  job  that  an  asymptotic  series  can  do  in  representing  the 
current  in  the  shadowed  parts  of  the  surface.  The  failure  of  the  L-K  asymptotic  scries  is  linked 
to  the  fact  that  the  current  in  the  shadow-zones  of  the  surface  cannot  ke  represented  by  an 
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asymptotic  scries  of  the  L-K  form.  To  prove  this,  recall  that  the  definition  of  an  asymptotic 
series  such  as  (2.11)  is  that  if  Sm  represents  the  partial  sum  of  the  first  m+  1  terms  *  then 


.  lim  C  f  ZV0,7)-sJ  «  0  (2.24) 

k0-*oo  L  J 

for  alt  values  of  m.  In  the  limit  as  kQ-+oo,  the  current  in  the  shadow-zones  of  the  surface  is  zero. 
Thus,  from  (2.24)  with  m=0,  S0  =  0  or  jjr )  =0.  For  m=  1,  (2.24)  yields 

lim  {k0L  )  =  lim  (^0Sj) 

kQ—foo  k0—*oe> 

but  =  j\lkg  because  j0  =0,  so  this  leads  to 

h  =  lim  (k0  L  )  (2.25) 

kQ-oo 

It  is  well  known  that  the  current  in  the  shadow  zone  has  the  form  of  a  creeping  wave  (11) 

L  -XCf  exp(  -j[k06  +  exp(  -jn^^kfl) 1/3  f*  K2/3^]}  (2.26) 

where  C(  are  the  launching  amplitudes,  6  is  the  distance  measured  along  the  surface,  pf  is  a  root 
of  the  Airy  integral  [11],  and  k  is  the  curvature  of  the  surface  at  the  distance  S  measured  along 
the  surface.  Thus,  in  the  limit  as  kQ-+oo,  each  term  in  the  scries  of  (2.26)  exhibits  an  exponential 
dependence  on  k0.  Furthermore,  since  the  real  part  of  this  dependence  leads  to  an  exponential 
decay  with  distance,  it  is  clear  that 


y'j  =  lim  (kQ  L  )  =  0 
k o 

and,  in  fact,  all  of  the  higher  jn 's  will  vanish.  Thus,  the  only  acceptable  asymptotic  series  for 
the  current  in  the  shadow  region  is  the  null  series.  The  reason  for  this  is  contained  in  the 
definition  of  the  asymptotic  scries,  e.g.  (2.24).  That  is.  the  only  scries  of  the  form 

m  — 

Sm  =  £  in(r  )lk0 
/!=() 

which  can  satisfy  (2.24)  for  all  m  in  the  shadowed  regions  of  the  surface  is  the  null  scries.  Note 
that  it  is  the  constraint  imposed  by  the  form  of  the  asymptotic  scries  which  dictates  the  end 
result.  Thus,  other  than  the  null  series  result,  the  current  in  the  shadow  region  does  not  have  an 

—  m  — 

*  Sm  =  I  J„(r 

n =t) 
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asymptotic  series  representation.  This  is  an  important  result  because  it  is  the  first  time  (to  the 
author  s  knowledge)  that  the  failure  of  an  L-K  series  in  the  shadow  region  has  been  both 
demonstrated  and  explained. 

While  the  L-K  scries  does  not  lead  to  an  exact  result  for  the  surface  current  density,  it  still 
holds  the  potential  for  providing  a  tractable  improvement  to  a  pure  geometrical  optics  solution. 
This  fact  is  demonstrated  by  Ansorge's  [12]  calculations  for  scattering  by  a  dielectric  sphere 
using  the  ray  optics  field  approach.  The  attractiveness  of  the  current  approach  as  developed 
here  is  due  in  large  part  to  the  fact  that  the  complete  L-K  representation  can  be  developed 
entirely  from  an  integral  of  known  functions,  i.e. 

/  10G{  |  Ar  | )  exp  (jkj .  A r  )dr,Q 

While  a  complete  L-K  series  development  for  this  integral  is  prohibitive,  it  should  be  possible 
to  obtain  the  terms  up  to  and  including  k~~  [9],  One  of  the  primary  advantages  of  obtaining 
the  k~l  and  k~ 2  corrections  to  the  k ®  asymptotic  expansion  of  this  integral  is  the  recovery  of 
some  of  the  cross  polarizing  properties  of  a  rough  surface  in  the  high  frequency  (but  not  optical) 
limit. 

Having  developed  the  Luneburg- Kline  expansion  for  the  current,  it  is  a  relatively 
straightforward  matter  to  find  the  scattered  field.  The  far-ficld  approximation  for  the  scattered 
magnetic  field  is 


—  -*  exp(  —jknR)  a  „  — 

HS{R)  =  -j  - — -  ksx  f  [k(TJ)exp{jks.r0)dr(o  (2.27) 

A 

where  R  is  the  distance  from  the  origin  on  the  surface  to  the  point  of  observation  and  ks 
specifiers  the  direction,  e.g.  R  =  Rky  Substituting  the  Luneburg-Kline  expansion  for  J  in  the 
above  and  assuming  an  interchange  of  the  summation  and  integration  yields 


Hs  =  -j 


exp(  -jk0R) 
4  itR 


A 

ksx 


V  1  r  i^rp) 

i  n  k 
rt=0  Ko 


exp [j(ks  -  kj)  .7y]  dr, 


The  integral  is  expanded  in  a  Luneburg-Kline  scries  as  follows; 


(2.28) 


I 


hiro) 


exp [j(ks  -  k() .  r0]dr,  =  I 
°  m=0 


nm 
i  in 


(2.29) 


where 
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hno  =,  U®  /- 
ko  00 


exp [/(*,  -  dr. 


(2.30a) 


A„1  =  lim  / 


expO(*j  -  ki)'~r0]  *t0  ~  hn 


(2.306) 


V  =,  i«m  ko  J  JLr~  exp0'(^  -  ki)-~ol  dnn  -  S 


m= 0  k„ 


(2.30c) 


so  that  //5  is  given  by 


W5=  -/ 


,  JkOR  A 

=— ——  fcj  X  X  X 

4ff/*  n=0  m-o 


The  geometrical  optics  field  is  proportional  to  hQO  which  is  determined  by  ja  only.  The  part  of 
Hs  which  varies  as  A:”1  depends  on  both  j0  and  yj,  the  part  that  varies  like  k~ 2  requires 
jQ,  j\,  and  y'2,  and  so  forth  for  higher  order  terms.  It  should  be  noted  that  with  the  physical 
optics  approximation,  one  obtains  only  those  contributions  from  jQ.  However,  it  is  also  known 
that  only  the  frequency  independent  term  is  always  accurate  and  the  above  series  shows  why  this 
is  true,  i.e.  the  y'j  contribution  to  the  /c~*  term  is  ignored,  the  y'j  and y'2  contributions  to  k~ 2  term 
are  ignored,  etc.  Inclusion  of  these  y'j,  y'2,  etc.,  terms  will  definitely  improve  the  frequency 
dependence  of  the  result  relative  to  the  physical  optics  approximation.  There  is  another 
advantage  to  including  these  terms.  Because  of  the  asymptotic  nature  of  the  series  being  dealt 
with  here,  each  term  by  itself  will  be  accurate  only  over  a  limited  range  of  observation 
directions.  For  example,  the  part  of  Hs  contributed  by  jQ  is  truly  only  accurate  about  the 
specular  direction,  i.e.,  ks  =  kf  and  ks ^  =  -kj  .  Thus,  if  more  jn  terms  are  included  it  should 
make  the  resulting  expression  for  Hs  more  accurate  in  directions  away  from  specular.  Clearly, 
the  calculation  of  these  jn  s  can  lead  to  definite  improvements  in  existing  models  and  is  therefore 
strongly  recommended. 
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